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Research  has  been  carried  out  in  the  areas  of  (1)  Saturation  spectro¬ 
scopy  including  effects  of  level  degeneracy,  (2)  Heating  and  cooling  of 
vapors  using  collisionally-aided  radiative  excitation,  (3)  Creation  of 
electronic  state  coherences  in  laser-assisted  collisions,  (4)  Two-level 
problem  plus  radiation  pulse,  (5)  Photon  as  catalyst  effect  and  (6) 
Collisional  processes  in  4-wave  mixing  experiments. 

(1)  Saturation  Spectroscopy  Including  Effects  of  Level  Degeneracy 
(C.  Feuillade,  P.  Berman). 

Owing  to  a  favorable  resonance  transition  frequency,  Na  has  been  the 
favorite  choice  of  experimentalists  in  laser  spectroscopic  studies.  The 
fine  and  hyperfine  structure  of  Ha  leads  to  a  multitude  of  levels,  even  in 
the  Na  ground  state.  There  have  been  no  rigorous  calculations  that  properly 
incorporate  the  effects  of  fine  and  hyperfine  structure,  collisional  effects 
and  optical  pumping  effects  with  Na  as  the  active  atom  in  a  laser  spectro¬ 
scopy  experiment.  However,  it  is  clear  that  optical  pumping  of  the  ground 
state,  in  particular,  can  severely  modify  the  laser  spectroscopic  line 
shapes. 

Due  to  the  fundamental  Importance  of  the  Na  system  in  laser  spectro¬ 
scopy,  we  have  begun  a  project  to  include  all  fine  and  hyperfine  structure 
of  the  3S,  3P  and  4D  levels  of  Na,  interacting  with  two  laser  fields. 

Both  steady  state  and  transient  solutions  will  be  sought,  to  clearly 
isolate  the  effects  of  optical  pumping.  Eventually,  collisional  effects 
will  be  Included. 

X 

The  first  stage  of  this  calculation  has  now  been  completed.  Using 
both  an  Irreducible  tensor  and  standard  (m-basis)  representation  for  atomic 
density  matrix  elements,  we  have  derived  expressions  for  the  probe  absorption 
line  shape  when  a  pump  field  of  arbitrary  strength  and  polarization  drives 
a  transition  between  two  levels  (each  containing  a  number  of  degenerate 
magnetic  sublevels)  and  a  probe  field  of  arbitrary  strength  and  polari¬ 
zation  drives  a  coupled  transition.  In  essence,  the  calculation  is  one 
describing  the  saturation  spectroscopy  of  three-level  systems  including 
effects  of  level  degeneracy.  The  probe  absorption  line  shape  is  calculated 
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for  each  velocity  subclass  of  atoms  but,  at  this  stage  o.  our  work,  no 
average  over  the  atomic  velocity  distribution  has  been  included.  Results 
are  presented  for  weak  probe  fields  as  a  function  of  pump  field  strength, 
polarization  and  detuning.  Effects  of  optical  pumping  are  Included. 

The  results  clearly  display  the  effects  of  pump  field  strength  (Rabl 

splittings) ,  pump  field  polarization  (position  and  number  of  Rabi  split 

resonance  peaks)  pump  field  detuning  (position  of  resonance  peaks)  and 

optical  pumping  (relative  strengths  of  the  resonance  peaks) .  We  have 

shown  that  it  is  advantageous  to  use  the  m-basis  rather  than  the  irreducible 

tensor  representation  provided  that  the  laser  fields  are  either  circularly 

or  linearly  polarized,  if  collisional  effects  are  unimportant.  We  have 

also  shown  how  to  predict  the  position  of  the  resonance  peaks  using  a 

2 

dressed-atom  approach. 

As  they  stand,  the  calculations  can  be  used  to  describe  the  inter¬ 
actions  of  laser  fields  with  an  atomic  beam.  The  next  step  in  the  calcu¬ 
lation  will  be  to  include  an  average  over  an  atomic  velocity  distribution 
so  that  laser  field-atomic  vapor  interactions  can  be  properly  modeled. 

(2)  Heating  and  Cooling  of  Vapors  Using  Collisionally-Aided  Radiative 
Excitation  (E.  Giacobino,  P.  Berman). 

3 

Several  years  ago  ,  we  predicted  that  cooling  or  heating  of  an  atomic, 
vapor  could  be  achieved  using  Collisionally-Aided  Radiative  Excitation 
(CARE).  In  Prof.  Stroke's  laboratory,  we  are  now  trying  to  carry  out  an 
experiment  of  this  type.  The  reaction  under  investigation  is 

Na  (3S1/2)  +  X  +  -W  +  Na  (3P1/2)  +  X 
where  X  is  a  rare  gas  atom. 

The  energy  defect  between  the  photon  energy  -fifi  and  the  3P^2  “  ^Sl/2 
transition  frequency  is  provided  by  a  correspond.ing  change  in  the  trans¬ 
lational  energy  of  the  Na  -  rare  gas  system.  To  probe  this  energy  change, 
the  velocity  distribution  of  the  excited  state  Na  atom  Is  monitored  using 
the  transition  to  the  AD  state.  Calculations  were  made  which  indicated 
that  heating  of  the  Na  should  be  detectable  by  this  scheme  using  a  positive 


energy  defect  and  heavy  rare  gas  perturbcrs.  Experimental  confirmation  of 

/«* 

the  heating  effect  has  now  been  achieved.  With  Ar  or  Xe  rare  gas  per- 
turbers,  the  excited  state  Na  velocity  distribution  was  considerably 
broader  than  the  thermal  one,  while  with  He  perturbers  [the  "light"  He 
takes  up  the  excess  energy  rather  than  the  "heavy"  Na]  the  normal  thermal 
width  was  observed. 

We  are  still  interested  in  producing  a  macroscopic  heating  or  cooling 

of  the  vapor  using  CARE.  It  may  be  possible  to  produce  such  an  effect 

using  a  high  density  sodium  cell.  In  that  case,  CARE  is  produced  via 

resonant  collisions  which  have  a  relatively  large  cross  section.  Order  of 

c  * 

magnitude  calculations  have  been  carried  out  which  indicate  that  measur¬ 
able  temperature  differences  can  be  achieved  with  Na  densities  of  order 
15  3 

10  /cm  .  An  experimental  test  of  the  predictions  is  envisioned. 

(3)  Creation  of  Electronic  State  Coherences  in  Laser-Assisted 
Collisions.  (P.  Berman,  E.  Giacobino). 

We  have  proposed  a  new  method  for  generating  electronic  state  co¬ 
herences  using  laser  assisted  collisions.  Both  Collisionally-Aided 
Radiative  Excitation  (CARE)  and  Radiatively-Aided  Inelastic  Collisions 
(RAIC)  can  be  used  to  generate  such  coherences.  Two  pulsed  laser  fields.  ' 
are  incident  on  two  atoms  undergoing  a  collision. 

In  CARE  (Fig.  la). 
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a 


V 


A' 


Fig.  la 

*  An  asterisk  on  a  reference  indicates  that  the  reference  is  appended  to 
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this  report. 


the  reaction  is 


A1  +  Al*  +  ^  **“  Ax(23)  +  Alt 

where  A  and  A'  are  the  atoms  undergoing  the  collision  and  the  notation 
A(iJ)  is  meant  to  indicate  a  coherence  between  states  i  and  j.  The  fre¬ 
quency  difference  -  ft  must  be  close  to  the  resonance  frequency  for 
the  coherence  to  be  generated. 

In  RAIC  (Fig.  lb), 


Fig.  lb 

a  typical  reaction  is  A£  +  A|,  +  -fift  +  A^  +  A  *  (2*3*)  producing  a 

2' 3'  coherence  in  atom  A*.  This  is  a  somewhat  novel  way  to  produce, 
electronic  state  coherences.  In  certain  cases,  it  may  be  possible  to 
produce  coherence  between  states  of  opposite  parity,  which  could  then 
radiate  sum  frequency  radiation.  This  process  is  under  investigation. 

(4)  Two-level  Atom  Plus  Radiation  Pulse  (E.  Robinson, 

P .  Berman) . 

Research  continues  in  the  fundamentally  important  problem  of  a  two- 
level  system  coupled  by  a  radiation  pulse.  In  the  large  detuning  limit, 
we  were  able  to  show  that  certain  classes  of  coupling  pulses  having  the 
same  asymptotic  Fourier  transforms  will  yield  transition  probabilities 
that  are  related  to  each  other  by  a  simple  scaling  transformation. 

Methods  for  evaluating  the  transition  probabilities  lit  the  large-detuning 
limit  have  also  been  developed,  and  we  are  trying  to  compare  our  results 
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with  those  of  other  authors.  '  In  particular,  if  we  would  like  to  under- 

O 

stand  the  details  of  Crother's  calculations  to  see  if  they  can  be  applied 
to  the  general  theory  of  laser-assisted  collisions. 

Work  continues  on  an  eigenvalue  method  for  solving  the  two-level 
problem.  Eigenvalue  expansions  for  the  probability  amplitudes  have  beer, 
obtained,  approximate  positions  of  the  eigenvalues  (i.e.  those  field 
strengths  leading  to  zero  transition  probability)  determined,  and  ex¬ 
pressions  for  the  transition  probability  to  third  order  in  the  detuning 
were  derived.*®*’*'* 

(5)  Photon  as  Catalyst  Effect  (A.  Lau). 

Using  experimentally  determined  values  for  energies,  transition 

moments  and  decay  rates,  the  rate  for  bound-continuum  transitions  in  I, 

12*  L 

resulting  from  the  photon  as  catalyst  effect  have  been  calculated. 

By  stimulated  emission,  the  laser  field  takes  an  bound  state  to  a 
virtual  intermediate  state  frota  which  a  transition  to  the  continuum  can 
occur  by  the  absorption  of  another  photon.  There  is  no  net  photon  absorption 
the  laser  field  acts  as  a  "catalyst". 

(C)  Collioional  Processes  in  A-Wave  Mixing  Experiments.  (P.  Berman) 

In  collaboration  with  J.  Lam  (Hughes  Research),  we  are  continuing  to 
attcuot  to  develop  a  consistent  theory  of  collision- induced  structure  in 

13 

4-wave  nixing  experiments  that  was  discussed  in  last  year's  Annual  Report. 

That  there  is  nuill  considerable  experimental  inte.rest  in  thic  area  is 

14 

evidenced  by  the  recent  work  in  Bloembergen' a  group. 

(7)  MiceeJ 3 eneous 

15  * 

Our  review  article  on  laser- esc is ted  collisions  has  been  published. 
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Abstract  The  physical  principles  underlying  the  com¬ 
bined  radiation  field  -  collisional  excitation  of  atoms 
are  reviewed.  A  discussion  of  both  collisionally- 
aided  radiative  excitation  ("optical  collisions")  and 
rndiatively-aided  inelastic  collisions  ("radiative 
collisions")  is  presented. 


INTRODUCTION 

The  purpose  of  this  paper  is  to  present  a  simple  discussion 
of  atomic  transitions  induced  by  the  simultaneous  action  of 
a  laser  field  end  a  collision. 

Consider  a  reaction  of  the  general  form 

Ai  +  Bi  +  ^  Af  *  Bf*  '  (1) 

where  A^  f  and  Bi  f  are  internal  states  of  two  atoms  A  and 
B  undergoing  a  collision  and  f]  is  the  frequency  of  an 
applied  radiation  field.  If,  in  the  absence  of  the  colli¬ 
sion,  one  finds 

aa  +  tin  ki 

Bi  +  tin  •*>  Bi 

while,  in  the  absence  of  the  external  field,  one  has 
A.  +  B,  A.  +  j.  . 
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the  reaction  (l)  is  of  a  type  that  requires  the  simultan¬ 
eous  presence  of  both  a  collisional  interaction  and  exter¬ 
nal  radiation  field  if  either  or  both  final  atomic  states 
are  to  differ  from  the  initial  ones .  One  may  then  speak  of 
"laser-assisted  collisions"  or  "collisionally-assisted 
light  absorption".  These  are  processes  vhich  have  been  the 
focus  of  a  large  number  of  experimental  and  theoreti¬ 

cal^”1*^  investigations  in  the  last  decade.  In  this  work, 
we  discuss  the  physical  principles  underlying  such  reac¬ 
tions;  more  detailed  theoretical  treatments  may  be  found  in 
the  literature. 

Reactions  of  the  form  (l)  may  be  further  classified 
into  two  categories.  The  first  of  these  we  refer  to  as 
Collisionally-Aided  Radiative  Excitation^  (CARE)  and  has 
been  designated  by  others  as  "optical  collisions".  The 
CARE  reaction  is  easy  to  visualize  (see  Fig.  l).  An  atom  A 


Figure  1.  A  schematic  representation  of  the  CARE  re¬ 
action  A^  +  B^+  -fin  ■*  A ^  +  B^.  A  laser  field  of  fre¬ 
quency  ft  is  incident  on  atom  A  and  c&n  drive  the  1-2 
transition  when  atom  A  undergoes  a  collision  with  a 
ground  state  perturber  B. 

is  irradiated  by  a  laser  field  whose  frequency  Q  is  close 


enough  to  that  of  „ 
approximation  to  l. 
tuned  from  exact  i, 
compared  to  the  nu 
tion,  but  small  cv 
With  such  a  large  i 
to  excite  atom  A  } 
goes  a  collision  v 
field,  the  probabi 
hanced.  The  ener, 
atomic  transition 
ponding  change  in 
atoms . 

The  second  c 
to  as  Radiativeiy 
has  been  designnt 
or  "LICET  -  Laser 

•3 

fer".  Atoms  A  ar, 
Bi  and,  as  a  cono¬ 
field  interaction 
Af  and  Bf.  The  p 

The  transit' 
assumed  to  be  hi, 
in  the  absence  o' 
the  photon  a 3  pr- 
transition  A^  +  ■ 
section  will  be 
to  be  resonant  v 
and  final  compon 
significant  exci 
tions,  with  the 


— ****\ 


W  r  » .  ■«,.  V- 


simultan- 
and  exter- 
nic  states 
hen  speak  of 
isisted 
ave  been  the 
theoreti- 
this  work, 
ueh  reac- 
be  found  in 


lassified 
or  to  as 
E)  and  has 
3". 19  The 
•  An  atom  A 


o  CARE  re- 
of  fre- 
the  1-2 
with  a 


}  is  close 


COMBINED  RADIATION  FIELD  -  COLLISIONAL  EXCITATION 

enough  to  that  of  an  atomic  transition  for  a  tvo-level 
approximation  to  be  valid.  The  field's  frequency  is  de¬ 
tuned  from  exact  resonance  by  an  amount  A  vhich  is  large 
compared  to  the  natural  and  Doppler  vidths  of  the  transi¬ 
tion,  but  8 mall  compared  to  the  thermal  energy  divided  by-fi. 
With  such  a  large  detuning,  the  probability  for  the  field 
to  excite  atom  A  is  negligibly  small.  However,  if  A  under¬ 
goes  a  collision  with  atom  B  while  interacting  with  the 
field,  the  probability  for  excitation  can  be  greatly  en¬ 
hanced.  The  energy  mismatch  -HA  between  the  photon  and 
atomic  transition  energies  is  compensated  for  by  a  corres¬ 
ponding  change  in  the  translational  energy  of  the  colliding 
atoms . 

The  second  class  of  reactions  of  the  type  (l)  we  refer 

to  as  Radiatively-Aided  Inelastic  Collisions  (RAIC)  and 

17 

has  been  designated  by  others  as  "radiative  collisions" 
or  "LICET  -  Laser  Induced  Colli si on al  Excitation  Trans- 
fer".  Atoms  A  and  B  are  prepared  in  initial  states  A^  and 
and,  as  a  consequency  of  the  combined  atom-atom  and  atom- 
field  interactions ,  they  emerge  In  some  new  final  states 
Af  and  Bf.  The  process  is  depicted  schematically  in  Fig.  2. 

The  transition  between  initial  end  final  states  is 
assumed  to  be  highly  improbable  or  energetically  forbidden 
in  the  absence  of  the  applied  field.  Thus,  one  can  view 
the  photon  as  providing  the  energy  to  assist  the  inelastic 
transition  A^  ♦  *►  A^,  +  B^..  In  general  the  RAIC  cross- 

section  will  be  largest  if  the  photon  frequency  is  chosen 
to  be  resonant  with  the  energy  difference  between  initial 
and  final  composite  atomic  states.  However,  as  in  CARE, 
significant  excitation  can  occur  under  off-resonance  condi¬ 
tions,  with  the  energy  mismatch  again  compensated  by  a 
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Figure  2.  Atoms  A  and  B  undergo  a  collision  in  the  pre¬ 
sence  of  radiation.  The  field  frequency  fi  is  approxi¬ 
mately  equal  to  a  transition  frequency  in  the  composite 
AB  system.  The  RAIC  reaction  is  of  the  form 
Ai  +  B.  +  -nil  Af  +  Ex¬ 
change  in  translational  energy. 

Before  examining  CARE  and  P.AIC  in  greater  detail,  it  is 
useful  to  review  the  problem  of  the  interaction  of  a  radia¬ 
tion  pulse  with  a  two-level  atonic  system. 

ATOM  +  PULSE 

In  this  section,  we  examine  the  interaction  of  a  two-level 
atom  with  a  radiation  pulse  whose  electric  field  of  polari¬ 
zation  e  may  be  represented  by 

E(t)  =  eEo(t)cos(fit) , 

The  smooth  pulse  envelope  function  ^^(t)  is  assumed  to 
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vanish  as  t  ±°°  and  to  vary  slowly-  in  comparison  vith 
cos(fit).  The  difference  -  0)| ,  where  u>  is  the  atomic 
transition  frequency,  is  taken  to  be  much  less  than  (u  +  ft) , 
allowing  one  to  neglect  the  "anti-rotating"  components  of 
the  field.  For  an  atom  which  is  in  its  lower  state  1  at 
t  =  -  00 ,  we  seek  the  probability  that  it  is  excited  to 
state  2  following  its  interaction  with  the  pulse.  Taking 
the  atom-field  interaction  to  be 

V(t)  =  -J4*E(t), 

where  jj  is  the  atomic  dipole  moment  operator,  one  may  use 
Schrodingcr's  equation  to  obtain  the  time  evolution  equa¬ 
tions  for  the  state  amplitudes.  In  the  interaction  re¬ 
presentation  and  with  the  neglect  of  the  anti-rotating 
components  of  the  field,  one  finds 

=  -ix(t)  elAt  a2  (2aj 

&2  =  -ix(t)  e-iAt  &1%  (2b! 

where  A  =  ft-ui  is  the  detuning,  x^)  =  pEo(t)/2ft  is  the 
coupling  parameter,  and  p  =  <l|^re|2>  =  p*.  The  frequency 
X(t)  is  sometimes  referred  to  as  the  Rabi  frequency. 

The  problem  is  conveniently  described  in  terms  of  the 
following  parameters:  (l)  the  puJce  duration  T,  (2)  the 
frequency  f  =  t ) /x( t )  which  determines  the  frequency  com¬ 
ponents  characterizing  the  pulse,  (3)  the  natural  lifetimes 
of  states  1  and  2  which  are  taken  to  be  much  longer  than  T, 
Justifying  the  omission  of  decay  terms  in  Eqs .  (2),  (M  the 
detuning  A,  and  (5)  the  Rabi  frequency  As  a  simpli¬ 

fication,  we  set  f  =  T"1,  which  is  a  good  approximation  for 
smooth  pulses . 
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If  the  detuning  and  envelope  function  are  such  that 
|A|T  »  1,  the  pulse  contains  negligibly  small  Fourier 
components  at  the  frequency  needed  to  compensate  for  the 
detuning.  In  this  limit,  the  pulse  is  said  to  be  adiabatic. 

That  is,  the  excitation  probability  following  the  passage 
of  the  pulse  is  vanishingly  small,  i_.e_. ,  proportional  to 
exp(-2|A|T)  for  typical  envelope  functions.  It  is  inter¬ 
esting  to  note  that  the  excitation  probability  remains 
exponentially  small  regardless  of  field  strength  x^)*  re“ 
fleeting  the  fact  that  the  Fourier  components  needed  to  j 

effect  the  excitation  are  essentially  absent.  As  the  field  . 
strength  x(*)  increases,  the  excitation  probability,  which  j 

is  proportional  to  A2  =  j/”  x(t)dt|2  for  A2  <<  1,  exhibits  j 

some  type  of  saturation  behavior  for  A  >  1.  Thus,  with-  ( 

out  some  additional  interaction,  an  adiabatic  pulse  cannot  j 

appreciably  excite  the  atom.  The  "additional  interaction"  j 

can  be  provided  by  a  collision. 

CARE'  I 

Assume  that  the  atom  undergoes  a  collision  with  a  perturber  J 
during  its  interaction  with  the  adiabatic  radiation  pulse.  i 

This  collision  occurs  on  a  time  scale  T  (typically  j 

—12  ^  I 

10"  sec  for  the  thermal  atoms  under  consideration  here) 
which  is  short  compared  to  T  (typically  10"^sec).  The  j 

perturber  can  be  considered  as  providing  an  effective  time-  j 

dependent  potential  which  modifies  the  energy  separation  of  • 
states  1  and  2  in  a  transient  manner.  If 'EY  (t)  is  the  I 

1  I 

collision-induced  modification  of  level  i's  energy,  then  j 

the  instantaneous  transition  frequency  is 

I 

w(t)  *  u  +  V^gft) , 

( 
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where 

VLS(t)  =  V2(t)  " 

It  is  implicitly  assumed  that  V^C t)  1  V2(t),  as  is  gener¬ 
ally  the  case  if  levels  1  and  2  belong  to  different  elec¬ 
tronic  configurations  The  collisio*  does  not  have  suffi¬ 
cient  energy  to  couple  levels  1  and  2  in  the  absence  of 
the  field  (see  Fig.  l). 

The  effect  of  the  collision-induced  transient  vari¬ 
ation  of  the  transition  frequency  is  to  introduce  appre¬ 
ciable  Fourier  components  into  the  excitation  mechanism  at 
frequencies  up  to  u)c  =  t'1  »  T"1.  These  added  Fourier 
components  lead  to  a  new  contribution  to  the  excitation 
probability  vhich  is  much  larger  than  the  exp(-2|A|T)  term 
associated  with  the  atom- adiabatic  pulse  interaction.  This 
"collisionally-assisted"  contribution  leads  to  a  CARE  re¬ 
action  of  the  form 


perturber 
n  pulse. 
i 

here) 

The 

Lve  time- 
•ation  of 
i  the 
.  then 


i 

| 

i 


A1  +  B1  +  'fin  ■*’  A2  +  V 

The  state  amplitudes  now  evolve  according  to 

*1  =  -ix(t)exp[iAt-  )df  )&2  (3a) 

42  *  -ix(t)expf.iAt+i/JvLS(t,)dt']a1,  (3b) 

subject  to  the  initial  conditions 

a1(-°°)  =  1,  a2(-»)  =  0.  (3c) 

In  order  to  discuss  CARE,  it  Is  useful  to  again  refer 

to  the  various  time  scales  in  the  problem.  The  collision 

duration,  Tc(b,vr)  -  b/v^,  where  b  is  the  intact  parameter 

and  vr  the  interatomic  speed  associated  with  a  collision, is 

an  important  time  parameter.  Although  t  (b,v  )  varies 

c  r 
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from  collision  to  collision,  ve  can  define  a  representative 

time  T  =  t  (b  ,v  )  in  which  v  is  the  average  interatomic 
c  c  o  r  r 

relative  speed  and  is  an  impact  parameter  chosen  to 
guarantee  that  T  is  "representative".  Generally  speaking, 
b0  will  be  that  impact  parameter  for  which  the  phase 

CO  w 

/  ^  V^gCb.v  ,t)dt  takes  on  a  value  of  order  unity;  a  typi- 
cal  value  for  bQ  is  10~  cm.  The  dimensionless  parameters 
which  enter  our  considerations  are  | A | T  which  turns  out  to 
be  unimportant,  { A | Tc  which  critically  categorises  the  de¬ 
tuning,  x(t)T  which  represents  the  strength  of  the  atom- 
field  interaction  before  and  after  the  collision,  and  yT 
vhich  represents  the  strength  of  the  atom-field  interaction 
during  the  collision.  The  field  strength  x(t)  is  approxi¬ 
mately  constant  during  a  collision  and  x  represents  some 
characteristic  value  of  |x(t)|  for  the  pulse.  As  noted 
above,  tc/T  «  1. 

Weak  Fields:  yT  «  1 

For  weak  fields,  the  excitation  probability  can  be  calcu¬ 
lated  from  Eqs .  (3)  using  first-order  perturbation  theory. 
The  results  depend  critically  on  the  value  of  | A | Tc . 

If  | A | xc  «  1,  the  only  change  in  state  amplitude  ag 
during  the  collision  arises  from  the  level-shifting  term. 
The  collision  acts  to  provide  a  sudden  change  in  the  phase 
4>  of  a2,  given  by  <|>(b,vr)  =  V^b.v^tjdt.  This  im¬ 
pulse  destroys  the  adiabatic  response  of  the  two-level 
system,  and  gives  a  final  state  amplitude 

a2  *  -il£  X(f)  e~iAt'dt'  ♦  e1*  /“  xU'  )e"iAt'dt*  ] 

c 

=  -  2ilx(tc)/A]e_iAtc  ei<^^sin(({i/2)  , 

where  tc  is  the  time  at  which  the  collision  occurs.  Set¬ 
ting  lx(tc)|  =  X»  one  obtains  the  excitation  probability 


la2(b'vr*w)r 

and  the  correspond 
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|a2(b,vr,«>)  |2=  l»(x/A)28in2[<j>(b,vr)/2]  (1») 

and  the  corresponding  CARE  cross  section 

Oc(vr)  =  2n/“|a2(b,vr,«)  l^db  (5a) 

*  2(x/A)2{7ib2).  (5b) 

The  result  ( U— 5 )  is  known  as  the  "impact  limit"  since  Tc  is 
smaller  than  any  other  time  scale  in  the  problem.  The  im¬ 
pact  cross-section  is  independent  of  the  sign  of  A  since 
the  Fourier  transform  of  the  collision  interaction  is  fist 
over  the  range  of  A  represented  by  |a|tc  «  1. 

The  impact  result  can  be  viewed  in  an  alternative 
manner.  If  we  were  to  suddenly  interrupt  the  atom-radi¬ 
ation  pulse  interaction  at  any  time  t  ,  we  would,  on 

2C 

average,  find  a  population  [x(tc)/A]  in  the  upper  state. 

The  CARE  cross-section  is  equal  to  the  product  of  this  ex¬ 
citation  probability  and  the  collision  cross  section 
i= 

If  !a|tc  >  1,  the  phase  induced  in  a2  during  the  colli¬ 
sion  by  the  detuning  is  not  negligible,  and  the  impact  re¬ 
sult  is  not  valid.  As  we  have  seen,  one  consequence  of  the 

collision  is  to  shorten  the  relevant  time  from  T  to  T  ,  so 

_!  c 

that  appreciable  Fourier  components  up  to  Tc  are  intro¬ 
duced.  If  this  were  all  that  occurred,  one  would  expect  a 
CARE  transition  probability  that  varied  as  exp[-2|A| 
Tc(b,vr)J.  However,  there  is  an  additional  effect,  whose 
origin  may  be  seen  in  Fig.  3,  which  modifies  this  result. 

In  drawing  the  energy  levels  in  Fig.  3,  we  have  chosen 
VLS(t)  <  0;  the  case  for  arbitrary  V^Ct)  may  be  treated  by 
an  obvious  generalization  of  the  method  given  below. 
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Figure  3.  Energy  levels  of  atom  A  during  a  collision. 
"Instantaneous  resonances"  [ft  =  w(t)]  can  occur  for 
detunings  A  <  0  only;  for  A  >  0,  the  collision  detunes 
the  atomic  transition  further  from  resonance. 

The  CARE  cross  section  is  a  strongly  asymmetric 
function  of  A  when  |aJtc  >  1.  For  a  given  A  <  0,  colli¬ 
sions  can  always  produce  u(t)  =  ft  for  short  times  during 
1*8 

the  collision;  i^.e:. ,  the  systems  become  instantaneously 
resonant  with  the  field.  Such  times  are  labeled  t^  and 
t2  in  Fig.  3.  The  phase  of  a2  varies  rapidly  owing  to 
the  factor  exp(-iAt),  except  at  t^  and  t2,  where  the 
oscillation  is  suppressed  by  the  factor  exp  [i 
/ ^ ( t ' ) dt '  ] •  The  major  contributions  to  the  excitation 
amplitude  are  provided  by  these  times  of  stationary  phase. 
The  corresponding  CARE  cross-section  varies  as  an  inverse 
power  lav  in  |A|,  instead  of  the  exponential  that  charac¬ 
terizes  other  regimes.  The  fact  that  the  points  of 
stationary  phase  provide  the  major  contributions  to  a2(<*>) 
is  linked  to  the  condition  |A|tc  >  1.  That  is,  the  (pulse 
♦  collision)  doeB  not  contain  the  Fourier  components  at 
A  to  appreciably  excite  the  atom;  in  this  case  the  in¬ 
stantaneous  resonances  become  a  critical  feature.  In  the 
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impact  limit,  the  system  of  (pulse  +  collision)  does 
have  appreciable  Fourier  coefficients  at  A  so  that  the 
presence  or  absence  of  instantaneous  resonances  does  not 
affect  the  excitation  amplitude. 

In  contrast  to  the  A  <  0  case,  for  A  >  0  the  colli¬ 
sion  pushes  the  levels  further  away  from  resonance  (see 
Fig.  3).  The  net  result  of  this  level  displacement  is 
that  the  nonresonant  side  of  the  CARE  cross  section  falls 

off  exponentially  as  a  fractional  power  of  |A|t  ,  even 

19  49  50 

after  one  averages  over  impact  parameter.  *  * 

Thus,  the  CARE  cross-section  exhibits  a  marked  asymmetry, 
with  an  inverse  power  lav  dependence  on  | A |  on  one  side, 
end  an  exponential  decay  on  the  other.  A  typical  pro¬ 
file  is  shown  in  Fig.  !». 


|  Figure  U.  CARE  cross  section  as  a  function  of  |a|tc 

i  in  the  weak  field  limit,  xTc  a  l.OxlO-1*.  This  cross 

section  is  drawn  for  a  level-shifting  term  which  varies 
as  R~^  (R  is  the  interatomic  separation)  and  a  value 
b  *  l.lxl0"7cm  (see  Ref.  37). 


P.  R.  BERMAN  AND  E.  J.  ROBINSON 


COMBINED  li¬ 


lt  should  be  noted  that  CARE  cross-section  in  the 
veak-field  regime  can  also  he  obtained  using  traditional 

pressure  broadening  theories  of  linear  absorption  or 

.  .  1*9—52 

emission. 

Strong  Fields:  yT  >  1 

As  long  as  x  <  |A|,  the  previous  perturbative  treatment  is 

2 

valid  and  the  CARE  cross-section  is  proportional  to  x  • 

If  both  xT  >  1  and  x  >  |A|,  the  perturbation  theory  fails, 
and  a  strong  field  theory  is  required.  Space  limitations 
preclude  a  detailed  description  of  such  a  theory,  vhich 
is  conveniently  developed  using  a  quantized- field  - 
dressed-aton  approach,  but  ve  cite  some  of  the  results. 

For  x  >  |A|  and  y*  <  1  (vhich  implies  | A | T c  <<  1), 
one  is  still  in  the  impact  domain  since  the  collision 
time  rc  i3  the  shortest  time  scale  in  the  problem.  If 
the  atom  -  radiation  pulse  interaction  is  interrupted  at 
some  arbitrary  time,  one  vould  find  an  upper  state  popu¬ 
lation  approximately  equal  to  1/2  since  the  field  iB 
sufficiently  strong  (xT  >  l)  to  lead  to  equal  populations, 
on  average,  in  levels  1  and  ?..  (This  factor  of  1/?  should 
be  compared  with  the  average  population  (xM)  found  in 
the  weak  field  case).  Thus,  in  this  limit,  the  CARE 

O 

cross-section  is  approximately  equal  to  nb^/2,  indepen¬ 
dent  of  both  A  and  y* 

For  X  >  |A|  and  yi  >  1,  an  impact  theory  can  no 
longer  be  used.  During  the  collision,  the  field  ia  strong 
enough  to  lead  to  rapid  oscillations  (so-called  Rabi 
oscillations)  in  the  state  amplitudes.  Since  x  >  |A|,  . 
these  Rabi  oscillations  provide  the  dominant  phase  vari¬ 
ation  for  the  state  amplitudes;  the  effective  detuning  in 
the  problem  becomes  x  instead  of  |Aj.  There  is  no  possi¬ 
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bility  of  "instantaneous  resonances"  hero;  consequently, 
the  excitation  probability  varies  as  exp[-2xtc(b,vr) ] . 
Just  as  in  the  weaX- field  result  for  the  A  >  0  cose,  the 
CARE  cross-section  obtained  after  averaging  over  impact 
parameter  falls  off  exponentially  as  a  fractional  power 
Of  XT  -19 


A  typical  RAIC  reaction  of  the  form 

A.  +  B.  *  difi  -*■  A.  ♦  B. 
i  i  f  f 

is  illustrated  in  Fig.  2.  In  going  from  the  initial  to 
final  Btate  in  the  composite  AB  system,  a  number  of  inter¬ 
mediate  states  may  play  a  role.  However,  by  Bumming  over 
these  states  and  neglecting  the  effect  of  small  vari ations 
in  nonresonant  energy  denominators,  one  may  reduce  the 
problem  to  that  for  a  two-level  system  coupled  by  an 
effective  operator  U(t)  which  is  proportional  to  the 
product  of  the  radiation  field  amplitude  and  the  colli- 
sional  interaction.  Explicitly  we  write 

U(t)  «  li (X/w)  Vc(t), 

where  w  is  some  representative  frequency  denominator 
(w  »  x)  &nd  Vc(t)  is  the  collisional  interaction.  Since 
U(t)  =  0  in  the  absence  of  a  collision,  the  RAIC  inter¬ 
action  occurs  during  the  collision  only.  Thus  the  pulse 
time  T  plays  no  role  at  all  in  RAIC  -  the  relevant  time 
scale  in  the  problem  is  the  collision  duration  ic. 

The  Initial  and  final  state  amplitudes  for  the  com¬ 
bined  AB  system  (see  Fig.  2)  obey  the  equations 
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=  -i(x/5)Vc(t)expliAt  -  i/Jv^t'Jdt' j  af  (6a) 

4f  *  -i(x/«)Vc(t)exp(-iAt  +  (6b) 

vhere 

A  =  Htt  -  (Ef  -  Et) 

and  VLS  (t)  is  the  same  collisions!,  energy  level  shift  en¬ 
countered  in  CARE.  An  additional  contribution  to  level 
shifts  resulting  from  the  AC  Stark  effect  will  not  be  dis¬ 
cussed  here,  but  can  be  included  by  redefining  the  energy 
levels  at  their  Stark-shifted  values. 

Weak  Fields 

A  perturbative  treatment  is  valid  provided  |a^(t)|  «  1, 
in  which  case 

|af(°°)  |  =  |(x/w)/“o0  Vc(t)exp[-iAt  +  i/jV^ft’ )dt' ]dt| . 

if  |a|tc  «  1,  the  amplitude  for  excitation  is  independent 
of  A.  For  (a|tc  >  1,  one  again  finds  an  asymmetric  line  aw¬ 
ing  to  the  effects  of  instantaneous  resonances  which  occur 

53 

for  one  sign  of  A  but  not  the  other. 

The  functional  forms  of  V  (t)  and  V._(t)  determine 

c  Do 

where  the  maximum  RAIC  cross  section  occurs  as  a  function 

of  A.  In  a  typical  situation,  the  time-dependence  of 

Vfi(t)  and  V^g(t)  is  roughly  similar  and  the  maximum  RAIC 

cross  section  occurs  for  A  =  0.  However,  the  RAIC  maximum 

may  occur  for  A  4  0  if  the  duration  associated  with  V£(t) 

is  much  smaller  than  that  associated  with  VTO(t),  as 

33^ 

might  be  the  case  in  RAIC  charge  transfer  ,  where,  for 
an  interatomic  separation  R(t),  Vc(t)  «  exp[-CR(t)), 
while  V^ft)  «  lR(t))”n.  Under  these  conditions,  the 
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collisional  coupling  is  significant  only  vhen  t)  1b 
producing  a  relatively  large  variation  in  the  atoms' 
energy  levels.  Since  the  effective  level  separation  of 
the  composite  AB  system  is  no  longer  Ef  -  during  the 
time  that  the  collisional  coupling  occurs,  it  is  not  sur¬ 
prising  that  the  maximum  RAIC  cross  section  can  be  dis¬ 
placed  from  L  =  0. 

Strong  Fields 

To  get  some  idea  of  strong  field  effects,  consider  Eqs. 
(6)  in  the  limit  that  A  =  0  and  vith  the  level-shifting 
term  set  equal  to  zero.  In  that  case,  the  upper  state 
probability  varies  as 

|a2(b,vr,®)  |2  «  oin2[<{>R(b,vr)  ], 

vhere 

Vb,vr*  “  Vc(b,vr,t)dt. 

2 

The  RAIC  cross-section  is  equal  to  7rbR  vhere  bR  is  an  Im¬ 
pact  parameter  for  which  is  of  order  unity.  For  a 
power  law  potential  V  (t)  «  [R(t))_  ,  n  >  3,4>R  “  x/B 

p  ^  -I  •*  * 

and  bp  varies  as  (x/v  Y  vith  a  =  (n  -  l)-  .  The  RAIC 
k  r  .  2 

cross-section,  which  is  proportional  to  x  for  veak 

fields,  varies  as  ^2/(n-l)  as  the  square  root  of  the 

Intensity  for  n=3)  in  the  strong  field  limit.  For  strong 
fields,  owing  to  the  fact  that  bR  a  xa»  large  impact 
parameter  collisions  only  are  important  and  VLS  plays  a 
minor  role  for  such  collisions .  The  line  width  is  de¬ 
termined  by  the  inverse  collision  time  Tc^  =  vy/bR  * 

the  RAIC  profile  narrows  with  increasing  field 

strength. 
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COMBINED 


The  ratio  bD/b  can  be  used  as  a  measure  of  the  field 

rt  O 

strength.  If  b,,  >  b  ,  one  is  in  the  strong  field  region 
n  O 

since  the  upper  state  amplitude  saturates  at  radii  where 

the  level-shifting  effect  is  unimportant.  On  the  other 

hand,  for  bR  «  bQ,  the  collisional  coupling  can  not 

overcome  the  effects  of  level-shifting  and  a  perturbative 

U 

treatment  is  valid.  Typically  ,  the  transition  from  weak 
to  strong  field  occurs  for  field  strengths  of  order 
lO^W/cm'v  The  strong  field  effects  in  RAIC  and  CARE  are 
fundamentally  different.  In  RAIC,  the  upper  state  proba¬ 
bility  is  truly  saturated  by  the  fi eld-collisionai  inter¬ 
action.  In  CARE,  on  the  other  hand,  the  upper  state 
probability  amplitude  is  always  small  if  xTc  >>  It  is 
the  rapid  Rabi  oscillations  that  lead  to  a  decreasing 
CARE  cross  section  with  increasing  X  when  xtc  >:>  I  am i 
X/N  >  1. 

CONCLUSION 

We  have  presented  explanations  of  the  physical  processes 
underlying  combined  radiation  field-collisional  exci¬ 
tation  of  atomic  systems.  Alternative  approaches  could 
involve  a  "dressed-atesn"  description  or  a  molecular-state 
basis  calculation.  For  a  meaningful  theoretical  des¬ 
cription  of  CARE  and  RAIC,  one  must  use  accurate  inter¬ 
atomic  potentials  and  average  all  results  over  the  spatial 
and  temporal  extent  of  the  laser  pulse.  It  may  be  noted, 
however,  that  experimental  investigations  of  CARE  and 
RAIC  have  revealed  many  of  the  qualitative  features  dis¬ 
cussed  above. 

This  work  is  supported  by  the  U.S.  Office  of  Naval 
Research.  The  content  of  this  paper  is  based,  in  part,  on 
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Theory  of  electronic  state  coherences  produced  in 
combined  laser  field  -  collisional  reactions 


A 
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A  method  for  producing  electronic  state  coherences 
using  either  CARE  (Colllsionally-Aided  Radiative  Excitation 
or  "optical  collision")  or  LICET  (Laser-Induced  Collisional 
Energy  Transfer  or  "radiative  collision")  is  proposed.  Two 
atoms,  A  and  A',  collide  in  the  presence  of  two  pulsed  laser 
fields  having  frequencies  ft  and  ft^.  It  is  shown  that,  by 
choosing  (ft^  t  ft)  such  that  an  energy  conserving  transition 
can  occur  in  the  composite  AA'  system,  one  can  create  an 
electronic  state  coherence  in  the  A  or  A*  atoms.  The  co¬ 
herence  can  be  produced  between  states  of  the  same  or  of 
opposite  parity;  if  it  is  between  states  of  opposite  parity, 
coherent  emission  at  frequency  (ft  ±  ft^)  is  generated. 
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1.  Introduction 


Over  the  past  ten  years  ,  there  has  been  increased  interest  in 
reactions  of  the  form 

At  +  -*  Af  +  ,  (I) 

in  which  two  atoms  or  molecules  (A  and  A')  collide  in  the  presence  of  a 
laser  field  and  undergo  a  transition  from  Initial  state  (A^A^,)  to  final 
state  (A^A^ , ) .  It  is  assumed  that  the  energy  difference  between  the  final 
and  initial  states  in  each  of  the  atoms  as  well  as  In  the  combined  system 
is  such  that  no  transition  occurs  in  the  absence  of  either  the  collision 
or  the  radiation  field.  Since  the  simultaneous  action  of  both  the  collision  and 
radiation  field  is  needed  to  produce  the  transition,  we  refer  to  reactions 
of  type  (1)  as  Combined  Fleld-Collislonal  Transitions  (CFCT) .  It  has  been 
common  practice  to  further  separate  the  CFCT  into  two  categories. 

2 

The  first  of  these,  Colllsionally-Aided  Radiative  Excitation  (CARE)  , 

3 

(also  referred  to  as  "optical  collisions"  ), characterizes  processes  in  which 
the  internal  state  of  one  of  the  reactants  remains  unchanged.  That  is, 

CARE  reactions  are  of  the  form 

A1  +  Al'  +  ^  "*  A2  +  Ai>  •  (2) 

subject  to  the  condition  that  direct  excitation  of  atom  A  by  the  radiation 
field  in  the  absence  of  the  collision  is  energetically  forbidden.  In  CARE 
the  collision  provides  translational  kinetic  energy  to  compensate  for  the 
mismatch  between  the  field  frequency  ft  and  the  1-2  transition  frequency  of 
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atom  A.  The  CARE  cross-section  falls  to  zero  for  atom-field  detunings 
which  correspond  to  energies  that  are  larger  than  those  available  in  the 
colliding  atoms'  center-of-mass  frame. 

The  second  type  of  CFCT,  Laser- Induced  Colllsional  Energy  Transfer 
4 

(LICET)  [also  referred  to  as  Radlatively-  Aided  Inelastic  Collisions 

2  5 

(RAIC)  and  "radiative  collisions"  ],  characterizes  processes  in  which 
both  reactants  change  their  internal  states  [see  Eq.  (l)].  The  photon 
energy  -fift  is  approximately  equal  to  the  difference  between  the  final  and 
initial  state  Internal  energies  of  the  atoms.  As  such,  the  radiation 
field  provides  the  energy  to  drive  what  would  normally  be  an  energetically 
forbidden  Inelastic  collision.^ 

In  principle,  CFCT's  can  be  produced  using  either  CW  or  pulsed  laser 
fields.  To  date,  LICET  has  been  seen  with  pulBed  excitation  only,  while 
CARE  has  been  observed  with  both  CW  and  pulsed  laser  fields.  In  this  work 
we  shall  consider  situations  Involving  pulsed  laser  fields  only.  The 
transition  from  initial  to  final  state  occurs  during  the  application  of 
the  pulsed  laser  fields;  we  calculate  the  final  state  density  matrix 
elements  characterizing  the  atoms  immediately  after  the  passage  of  the 
radiation  pulses. 

The  calculation  of  CARE  and  LICET  cross  sections  with  pulsed  laser 
excitation  has  been  the  subject  of  a  large  number  of  theoretical  and  ex¬ 
perimental  studies.^  In  a  typical  calculation  or  experiment,  one  determines 
the  excitation  cross  section  as  a  function  of  the  strength  and  frequency 
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of  Che  radiation  field,  and  attempts  to  gain  some  information  concerning 
the  initial  and  final  state  A-A'  interatomic  potentials.  While  most  of 
the  calculations  and  experiments  have  concentrated  on  final  state  popu¬ 
lations,  there  has  been  some  recent  work  devoted  to  the  investigation  of 

the  possibility  of  creating  coherences  using  CFCT  with  pulsed  laser 

7  8 

fields.  It  has  been  shown  '  that,  in  cases  where  the  final  states  con¬ 
sist  of  a  number  of  degenerate  magnetic  sublevels,  CFCT  can  create  co¬ 
herences  among  these  sublevels.  The  existence  of  such  CFCT-lnduced 

magnetic  state  coherences  has  been  experimentally  established  for  both 

9  10 

CARE  and  LICET. 

In  this  article,  we  present  a  theory  of  a  new  type  of  CFCT,  in  which 
it  should  be  possible  to  produce  electronic-state  coherences  using  pulsed 
laser  excitation.  The  types  of  reactions  to  be  considered  are  conveniently 
represented  in  Figs.  1-3.  In  all  cases,  atoms  A  and  A'  are  undergoing 
a  collision  at  some  time  during  the  application  of  one  or  more  laser 
pulses.  The  laser  pulse  durations  T  are  much  greater  than  the  duration  of 
a  collision,  but  much  smaller  than  the  lifetimes  of  the  relevant  atomic  levels. 
Following  the  pulses,  we  wish  to  determine  whether  or  not  any  electronic 
state  coherences  have  been  created  by  the  CFCT.  By  convention,  a  state 
written  as  A(fg),  Implies  that  the  A  atoms  possess  a  coherence  between 
states  f  and  g. 

Corresponding  to  each  figure,  one  can  write  a  reaction  that  could 
possibly  lead  to  electronic-state  coherences.  For  Fig.  1,  the  CARE  re- 
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action  is 


A1  +  A',  +  ISft  ♦  A(12)  ♦  A^,  , 


(3) 


for  Fig.  2,  the  CARE  reaction  is 

Ai  +  Ai'  +  *0  ♦  UOj  ♦  A(2 ' 3 * )  +  Aj,  ,  (4) 

and,  for  Fig.  3,  the  LICET  reaction  la 

A2  +  A|,+  «)  +  ♦  Ax  +  A’  (2'3')  .  (5) 

A  calculation  of  the  actual  macroscopic  coherence.  If  any,  generated 
by  reactions  (3-5)  is  given  In  Secs.  Ill  and  IV,  following  a  discussion 
In  Sec.  II,  of  the  approximations  and  assumptions  of  the  theory.  In  Sec.  V 
we  indicate  methods  for  detecting  the  CFCT- Induced  coherence.  The  choice 
of  detection  scheme  depends  on  whether  the  coherence  is  created  between 
states  of  the  same  or  of  opposite  parity. 

The  level  scheme  of  Fig.  2  has  been  analyzed  for  CV  laser  fields 

11  12  13 

by  Bloembergen  and  coworkers  ,  Grynberg  ,  and  others  ;  in  some  sense, 

our  analysis  for  this  case  represents  the  pulsed  field  analogue  of  their 

work.  It  may  also  be  noted  that  our  theory,  while  containing  some  features 

found  in  calculations  involving  "adiabatic  following"  in  two  and  three- 
14 

levels  systems  ,  differs  considerably  in  spirit  and  content  from  those 
calculations.  There  is  a  somewhat  closer  connection  between  the  under¬ 
lying  theory  in  our  work  and  that  exposed  in  the  recent  article  of 
Agarwal  and  Cooper^,  but  the  approach  and  emphasis  of  the  two  calculations 
differ  appreciably. 
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11.  Notation,  Approximations  and  Assumptions 

The  stain  goal  of  this  article  is  to  illustrate  the  manner  in 
which  CFCT  using  pulsed  laser  fields  can  create  electronic  state  co¬ 
herences.  Consequently,  we  shall  make  a  number  of  assumptions  and 
approximations  to  simplify  the  theoretical  development.  Most  of  these 
assumptions  can  be  relaxed  in  a  more  complete  theory,  but  the  basic 
physical  content  of  the  theory  would  remain  unchanged. 

Before  listing  our  assusiptlons  and  approximations,  it  is  useful  to 
Introduce  some  notation.  The  atoms  A  and  A'  of  Figs.  2  and  3  are  sub¬ 
jected  to  two  laser  pulses  having  amplitudes  £ (t)  and  £j(t),  frequencies 
ft  and  ft^,  and  propagation  vectors  k  and  k^,  respectively.  (The  atoms  of 
Fig.  1  are  subjected  to  one  pulse  only).  The  laser  pulses  have  duration  T 
while  the  collision  duration  is  T£.  The  detunings  A  and  A^  refer  to 
atom-field  detunings  for  single-photon  transitions  [e.g.,  in  Fig.  2 
A  »  ft  -  u>21,  ^  •  ftj-  u>31;  and  x  and  refer  to  the 

Rabi  frequencies  [evaluated  using  values  of  |£(t)  |  or  |  £^(t) |  averaged 
over  the  pulse  duration]  for  single-photon  transitions.  Spontaneous 
decay  rates  are  denoted  by  y  and  collision  rates  by  T.  The  most  probable 
atomic  speed  is  u. 

The  following  assumptions  or  approximations  are  made:  (1)  The  laser 
pulse  duration  is  short  enough  to  neglect  any  spontaneous  emission  during 
the  pulse, 

YT  «  1. 
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(6) 


(2)  The  detunings  |A  I  and  jA^|  are  sufficiently  large  such  that 

M  T  »  1  ;  |Aj  T  »  1.  (7) 

This  condition  will  ensure  that  there  la  negligible  excitation  In  the 
absence  of  the  collision  provided  that  we  also  require  | A (  and  |A^|  to 
be  larger  than  the  Doppler  widths  associated  with  the  transitions,  i.e. 

|A|  »  ku  ;  | |  »  k^u  .  (8) 

(3)  The  collision  occurs  on  a  time  scale  short  compared  with  the  pulse 
duration 

T  /T  «  1  .  (9) 

c 

(4)  The  collision  aay  be  treated  In  the  impact  approximation.  Implying 
that 


r  tc  «  i  , 

(10) 

) A |  Tc  «  1  ;  |A2|  Tc  «  1  , 

(11) 

X  Tc  «  1  ;  X1  Tc  «  1  . 

(12) 

Conditions  (11)  and  (12)  simplify  the  mathematical  development,  but  are 
not  essential  to  the  theory.  (5)  The  ground-state  lifetimes  are  Infinite. 
(6)  Two-photon  processes  are  resonant  or  nearly-resonant.  For  Fig.  2 
this  assumption  takes  the  form 

\U  -  -  u>32|t  «  1  ,  (13) 

while  for  Fig.  3a  or  Fig.  3b,  It  is 

^  "  h)3'2'|T  <<  (14) 

where  the  sign  refers  to  Fig.  3a  and  the  "+"  sign  to  Fig.  3b. 


Conditions  (13)  and  (14)  ensure  tnat  tne  conerences  vnich  are  created 
do  not  "wash  out"  during  the  laser  pulse  duration  T.  In  parallel  with 

conditions  (13)  and  (14),  ve  must  require  that  the  Donoler  shifts  app¬ 
elated  with  two-photon  transitions  that  are  accumulated  during  the  laser 
pulse  be  negligible.  For  the  case  of  Figs.  2  and  3a,  we  assume 

|k  -  kJuT  «  1  (15) 

while  for  Fig.  (3b),  the  condition  is 

|£  +  SjuT  «  1  (16) 

(7)  In  order  to  be  able  to  calculate  the  coherences  using  perturbation 
theory,  we  assume  that 


[dxl2  *  IxJ2  ♦  IxXjh/^jn  «  i 

(17) 

“d12  Klxl  +  IXjI2)M|t«  i  . 

(8)  Finally,  we  must  require  that 

(18) 

TT  «  1 

(19) 

to  ensure  that  the  coherence  does  not  decay  during  the 

Together  with  Eq.  (7),  Eq.  (19)  Implies  that 

radiation  pulse. 

r  «  |a|  ;  r  «  |a1 i  . 


(20) 


III.  Generation  of  Coherences  Using  CARE. 

A.  Figure  1 

We  first  turn  our  attention  to  the  situation  depicted  in  Fig.  1. 
A  laser  pulse  of  the  fora 


E(R,t)  -  -r  £(R,t) 


A  C  -  fit) 


+  c.c. 


is  incident  on  a  vapor  containing  A  and  A'  atoms.  In  a  field  interaction 
representation  to  be  defined  below  and  in  the  absence  of  collisions, 
density  matrix  elements  of  atom  A  evolve  according  to 

P12  “  *  (2  Y2  +  1A)P12  “  iX  (t)(pll  "  p22^  (22a) 

P,,  ■  “  Y?P»  +  i[x(t)p19  -  X  (t)p9.  ]  o?v,\ 


(22b) 


where  y^  is  the  spontaneous  decay  rate  of  level  1, 
*  _  i (k*R-ftt) 


P12  ’  P12  6 

X(t)  -  <  2|px|l  >  £(R,t)/2fi, 


A  -  (2  -  w21  (25) 

and  px  is  the  x  component  of  the  atomic  dipole  moment  operator. 

In  writing  Eqs.  (22),  y^  has  been  set  equal  ro  zero  and  the  Doppler  shift 
has  been  neglected  relative  to  |A(  [Eq.  (6)].  In  the  impact  approximation 
[ iqs .  (10-12) J ,  collisions  are  incorporated  into  Eq. (22)  by  the  addition 
of  a  term  (-I^p^)  t0  the  right  l'and  side  Eq*  (22a)  fT^  is  the  rate 
at  which  collisions  destrov  ij  coherence***] .  If,  in  addition,  we  use 
condition  (17)  to  set  p^^  =  1  and  P?2-0  in  Eq.  (22a)  (perturbation  limit), 
the  appropriate  equations  to  be  solved  during  the  laser  pulse  are 
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®,  one  finds 


(30) 


This  result  is  consistent  with  conservation  of  energy  considerations.  If 
there  are  no  collisions,  then  the  excitation  probability  for  a  detuning  A 
is  determined  by  the  Fourier  transform  of  the  pulse  envelope  function 
evaluated  at  A.  If  |A|T»1  as  assumed  in  this  work,  the  field  does  not 
possess  the  proper  frequency  components  to  excite  the  atom  [typically,  the 
excitation  probability  varies  as  exp(-2|A|T)  «  l],  Once  collisions  are 
present,  however,  the  effective  time  scale  of  the  problem  is  changed  from  the 
pulse  duration  T  to  the  collision  duration  Tc«  The  collision  does  possess 
the  Fourier  components  to  excite  the  transition  (t  ^  >>  j A | )  and,  moreover, 
provides  a  mechanism  (change  in  center-of-mass  translational  energy  of  the 
colliding  atoms)  allowing  for  conservation  of  energy  in  the  total  (absorption 
+  collision)  process.  The  corresponding  collision  enhanced  excitation  pro¬ 
bability  is  indicated  in  Eq.  (30). 


It  might  be  thought  that  collisions  could  also  lead  to  a  non-negli- 


glble  macroscopic  coherence  following  the  radiation  pulse.  Collisions 
do,  in  fact,  produce  non-negliglble  coherence  each  time  they  occur  during 
the  radiation  pulse,  but,  on  averaging  over  the  various  times  during 
the  radiation  pulse  at  which  they  occur,  one  finds  that  the  coherence  is 
negligible.  Formally,  the  result  follows  directly  from  Eq .  (26a),  l.e. 


Pl2(T+)  -  -ie‘1AT 


.+  ,T+  -(f  Y2+r12)(T+-t’) 


eiAt  Y(t’)dt' 


(32) 


Using  Eqs.  (6)  and  (20)»  and  letting  t+  ■+  “  in  the  Integral,  we  obtain 
a  coherence 


*  /-T+1  _  4  -1AT 

Pl2 (T  )  •  -ie 


T  x* 


(t  * )eiAt  dt ' 


which,  for  | A  j  T  >>  1,  is  negligibly  small. 


(33) 


B.  Figure  2 

In  order  to  produce  coherences  which  do  not  dephase  during  the 
radiation  pulse,  one  can  use  the  level  scheme  depicted  in  Fig.  2.  Levels 
2  and  3  are  assumed  to  have  the  same  parity  which  is  opposite  to  that  of 
level  1.  Laser  fields  of  the  form 


E.R,t)  - 

V  * 


2. 


6  (.R ,  t )  * 


i£-  ft  -n  t ) 

*■  c .  c 


i  €,(  *  e 


-n,t) 


c .  t 


(34a) 

(34b) 


are  simultaneously  incident  on  the  vapor.  We  wish  to  calculate  the  co¬ 
herence  p,-(T+)  using  perturbation  theory.  In  a  field  interaction  re- 


representation  defined  by 


i(k*R-flt) 


p12  *  p12  6 

(35a) 

r  .  p  .K^-R-^t) 

p13  p13 

(35b) 

pi,  «  p„ 

^23  y23 

(35c) 

density  matrix  elements  evolve  according  to 

-  ;  a  \  ~  ,/(t)lc  -«  ,)  t  t  (36a) 


-  L  i  ^  ~  r"  *■  L  A<v) 


V  (36b) 


z  ^ 


,  .  [  i  (1W,)  <■  l’.» 


c  -  If  \ 

•  Jv  vS  .  J  / 


where 


(36d) 


Av  -  -k*v .  A„  »  At-u,,  -k,*v 


and 


xz;,^ 


21  v  •  v  “1  31  1 

<£  i\  p,  M  >  z  * 

lt>  £,(£0/** 


X3,  t)  -  <  M  p 


(37) 


(38) 


Equations  for  Pj2  and  need  not  be  written  since  populations  in  these 
states  created  by  CARE  excitation  will  not  appreciably  affect  when 
condition  (18)  holds.  To  calculate  P2j,  we  set  P22  *  ^33  "  p32  "  p23  " 
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and  ■  1  in  Eqs.  (36a)  and  (36b),  owing  to  the  perturbation  limit. 
With  these  simplifications,  Eqs.  (36)  reduce  to 


^3 


—  a.) 

A/ 

<?,*  ~ 

l  x*,(t) 

(39a) 

-  C  * 

■  rl%  *  a’ 

•  L  X-3.lt) 

(39b) 

-  t  iO'1 

♦Y\)  rf1% 

«-  HA  l 

-Av)]  <?i4  V 

l  xz>t)^i3  -  t 

(39c) 

(39d) 

The  calculation  of  the  coherence  p2^  parallels  that  for  p 22  given 
above  in  Sec.  A.  If  there  are  no  collisions,  the  coherence  created  is 
negligible.  With  collisions  present,  Eqs.  (39a)  and  (39b)  are  Integrated 
by  parts  twice  and  the  resulting  values  for  p 12  and  (neglecting  terms 

o  o 

of  order  x/|A|  or  x’/|A'|  )  are  substituted  into  Eq.  (39c)  which  is  then 
Integrated  to  yield 

^  {  L  Uy>vY»  )  *'r,e  ♦  } 


q2^-r)  * 


-  z 


C  (-  l  t  V-  '"S3  r  ^  (  r+’-  C  V. 


x  <L*  c 


(40) 


In  arriving  at  Eq.  (40) «  ve  have  used  the  fact  that 


A1  -  A  { 
v  v\  «  1 


which  follows  from  Eqs.  (7)>  (13)  and  (15).  The  last  two  terms  in  Eq. 
(40)  may  be  integrated  by  parts  to  give 

*  i 

'  [.  t  (W  *3  )  *-  +  )  ]  J  ^  ^  ^ 


Combining  Eqs.  (42)  and  (40),  using  conditions  (6),  (13)  and  (15)  to 
set  the  exponential  appearing  in  the  integrand  equal  to  unity  and  taking 
T+  -►  one  finds 


-7  ^  .-1  n 

l  ,i  *  »  »*  '  1  z* 


i— ’  (  C*.0  Jt  (43) 


where  the  spatial  dependence  of  the  fields  has  been  made  explicit.  Equation 
(43)  gives  the  coherence  created  by  CARE. 
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! 


The  result  (43)  can  be  viewed  as  the  pulse  analogue  of  pressure- 

11  12 

induced  extra  resonances  discussed  by  Bloembergen  et  al  ,  Grynberg 
13  «/ 

and  others.  A  steady-state  value  for  calculated  assuming  a 

repetition  rate  for  the  radiation  pulses  has  the  same  structure  as  one 

of  the  terms  contributing  to  given  in  the  works  cited  above.  ~  ' 

Inere  is  an  additional  term,  however,  which  appears  in  the  CW  calculation 

that  does  not  appear  in  our  pulsed  version.  This  additional  term  does 

not  vanish  in  the  absence  of  collisions.  Such  terms  appear  typically 

in  steady-state  theories  (for  example,  a  two-level  atom  driven  by  an 

2  2 

off-resonance  CW  field  acquires  a  steady-state  population  x  /A  for  any 
A  ).  However,  even  though  such  terms  appear,  one  must  probe  the  aystem 
on  a  time  scale  shorter  than  1/|A[  in  order  to  Isolate  their  contribution, 
(i.e.  to  "see"  the  population  in  the  two-level  example  mentioned  above 
one  must  turn  off  the  field  in  a  time  which  is  Bmall  compared  with 

N'1). 

It  might  also  be  noted  that  the  coherence  p^  can  be  seen  in  the 
CW  experiments  only  if  |Av~A^|<  (-j  y^  +  ^3).  The  pulsed  version 
creates  coherences  over  the  much  larger  range  of  detunings  |A-A^|  <  T  \ 
Methods  for  detecting  the  coherence  are  discussed  in  Sec.  V. 


IV.  Creation  of  Coherences  Using  LICET  (Figure  3) 

In  Fig.  3a  or  Fig.  3b,  a  coherence  *8  generated  by  the 

LICET  reaction 

f\  A  +■  "fe  n  +  "h  -TL  <  ^  A ,  r  A  (  2.  3  )  m  (44 

Once  created,  the  coherence  oscillates  at  frequency  s  ft 

(Fig.  3a)  or  2 ,  -  ft^  +  ft  (Fig.  3b)  so  that  it  is  possible  to  generate 
coherences  at  either  the  sun  or  difference  frequencies  of  the  applied 
laser  fields.  (In  CARE,  it  is  also  possible  to  generate  a  coherence  at 
the  sum  frequency  if  one  starts  with  atom  A  in  an  excited  state). 

To  obtain  equations  describing  the  time  evolution  of  atomic  state 
density  matrix  elements  during  the  passage  of  the  laser  pulses  which  pro¬ 
duce  the  LICET  reaction,  it  is  useful  to  note  an  important  difference  be¬ 
tween  CARE  and  LICET.  In  CARE,  non-negligible  excited  state  density  matrix 
elements  oi  atom  A  (see  Fig.  2)  are  produced  during  the  laser  pulses  even 
in  the  absence  of  collisions.  However,  in  the  absence  of  collisions,  these 
excited  state  density  matrix  elements  adlabatically  follow  the  field  and 
vanish  at  time  T+  immediately  after  the  passage  of  the  laser  pulses.  Colli¬ 
sions  break  this  adiabatic  following  and  lead  to  non-vanishing  excited  state 
density  matrix  elements  at  time  T+.  Consequently,  the  equations  that  de¬ 
termine  the  time  evolution  of  excited  state  density  matrix  elements  in  CARE 
contain  a  collision-independent  contribution  (the  adiabatic  -  following  term) 
plus  a  colllsional  contribution.  In  general,  the  collislonal  term  can  depend 
on  the  laser  field  strength  and  the  atom-field  detuning;  however,  in  the 
Impact  approximation,  the  collislonal  contribution  is  independent  of  the 


field  variables  and  can  be  represented  simply  by  the  rates  T 
Eqs.  (36). 


found  in 


On  the  other  hand,  atom  A'  excited-state  density  matrix  elements 
produced  in  LICET  (Fig.  3)  result  solely  from  the  combined  radiation  field- 
collisional  excitation;  they  are  neglibibly  small  during  the  passage  of 
the  laser  pulses  if  no  collision  occurs.  The  coherence  ^ i  can  be  created 
directly  during  a  collision;  the  average  value  for  P2t2»(R,T  )  produced  by 
this  "direct”  excitation  channel  is  calculated  below.  The  coherence 
P2,2i(R>T  )  can  be  produced  by  an  "indirect"  process  involving  (a) 

the  generation  of  populations  P2<2'  or  P3131  by  LICET  at  some  time  during 
the  laser  pulses  followed  by  (b)  the  laser  fields,  acting  in  the  absence 
of  collisions,  to  generate  the  coherence  #2* 3'  *ro°  the  difference 
(Oyy  -  P2f2i)«  In  this  Indirect  process,  the  population  and  coherence 
are  created  sequentially  during  the  same  laser  pulses.  The  contribution 

r<J 

to  P0,v  from  the  indirect  process,  while  easily  calculable,  can  be 
neglected  relative  to  that  of  the  direct  process  when  condition 
(18)  holds. 

In  light  of  the  above  discussion,  we  proceed  to  obtain  the  LICET-ln- 

duced  coherence  p_,,,  as  follows:  (1)  the  value  of  p- ,,,(b,v  ,t  ,R)  re- 
11  1  J  r  c 

suiting  from  a  single  collision  characterized  by  impact  parameter  b  and 
relative  speed  v^  is  calculated  for  a  collision  occurring  at  a  time.  t£ 
and  position  R.  (2)  From  this  value  for  P2 1 31 (b,vf,tc,R) ,  we  find  the 
average  2' -3'  coherence  following  the  laser  pulses  as 
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f  ( ai  *)  -  rvA  ^  an  b  ^  ^ 


X  5  Jtc  ?2.3.  (  fe,vk  ,ttJR) 

T  - 

(45) 

where  N  is  the  density  of  A-atoms  and  W(v  )  la  the  relative  speed  distri- 
a  r 

butlon. 

mJ  4 

The  calculation  of  p£ ' 3' ^»vr*tc»R)  Ik  most  conveniently  carried  out 
using  state  amplitudes  rather  than  density  matrix  elements.  The  laser 
fields  Incident  on  the  atoms  are  again  given  by  Eq.  (34).  When  a  collision 
occurs  during  the  laser  pulses,  the  state  amplitudes  change  as  a  result  of 
three  effects. ^  First,  there  Is  a  shifting  of  the  levels  owing  to  the 
light-shift  operator;  for  the  perturbation  calculation  of  this  paper,  this 
term  can  be  neglected.  Second,  there  is  a  shifting  of  the  levels  owing 
to  a  collislonal  operator.  The  collisional  level  shifts  become  important 
for  collisions  with  impact  parameters  less  than  or  of  the  order  of  some 
critical  Impact  parameter  bQ  (typically  of  the  order  of  the  Welsskopf 
radius  associated  with  theories  of  pressure  broadening).  For  collisions 
with  b  >  bo>  the  collislonal  level  shifts  can  be  ignored;  for  collisions 
with  b  <  bQ  there  is  a  rapid  phase  variation  of  Pj.^i  (b.v^t^.R)  with  b, 
leading  to  destructive  interference  in  the  integral  (45).  Thus,  in  calcu- 
lating  p2ij,(R,T  ),  we  should  aet  the  lower  limit  of  the  b  integration 
to  bQ  and  neglect  the  effects  of  the  collisional  shift  operator  in  treating 
collisions  with  b  >  bQ.  Third,  there  is  the  LICET  transition  operator 
which  couples  the  initial  and  final  states  of  the  AA*  system  and  gives 


rise  to  the  coherence 

Matrix  elements  of  the  LICET  transition  operator  between  the 
initial  state  |2  1*  >  and  final  state  <  lf'|  (f'  -  2'  or  3’)  are  given 


V  U.  V,  ,tc  ,  R  ;  t  )  --  0/2*] 


<  1  f "(  (  ee'Xee'l  U  L  |  2  1  '> 

*  (^ax  *  ^e'lO 

<w'l  ti  ee'Xee'l  f>;  |  2l'> 

^  C  +  “V* '  ) 


(46) 

_  /  / 

•T  =  ? 

(47) 

ll  is  the  AA’  collisional  Hamiltonian  evaluated  along  the  classical 
collision  trajectory  R^f(b,vi>t),  and  the  sum  is  over  all  intermediate 
states  e  and  e*  [primed  (unprimed)  quantities  refer  to  atom  A  (A')]. 

Note  that  T^,  represents  the  combined  action  of  the  radiation  field 

and  collision  and  vanishes  if  either  &  ft  °r  2C  is  *ero.  In  terms  of  the 

matrix  element  (46),  the  final  state  amplitude  a^ , (b,vr,tc,R) ,  calculated 
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in  lowest  order  perturbation  theory  and  in  the  impact  approximation 
[Eqs.  (10-12) J,  is  given  by 


al5.  U.M*.  ,tt^) 


V( 


b  v , 


where 
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A,  if) 


\  T 

A  -  k  •  V  -  XL  -  (^2-,'  _u,ii  )  -W-V 

-7  -7 

A(-k,'V  -  -  (U»v,-  -‘O-.I  )  -  » 


(48) 


(49) 


Lb  vA,lt,S)  -  5  ;  Tt,L,v-.L,‘i  °dt- 


(50) 


and  t~  and  t+  represent  times  before  and  after  a  collision  centered  at 
c  c 

fJ 

t  •  t  .  The  value  of  p, produced  by  this  collision  is 
c  *  ^ 


9  e 


_  a  Avci'  ntt 


Ta,  v  (b  ,  t)  , 


i  1  i 


(51) 


where 


t  .  b v  -  r.,  n.nr  . 


(52) 


In  carrying  out  the  Integration  (45),  we  can  take  the  exnonential  factor  ir. 

Eq.  (51)  to  be  constant,  owing  to  conditions  (14)  and  (15).  Settlne  t  *  n 

c 

without  loss  of  generality,  we  find  the  value  of  ^,.,(R,T+)  following  the 

4-  J 
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pulses  fro*  End.  (45)  and  (51)  to  he 


=>= 

q  .  (  jC  T  M  -  ^  1  ]  W  ( VA )  cl 

»  iTT.dt,  VUA,'..?) 

(53) 

where  T^,^,  Is  calculated  from  Eqs.  (52),  (50)  and  (46). 

Without  going  Into  a  detailed  evaluation  of  Eq.  (51),  one  can  note 
several  general  features  of  the  result.  First,  the  magnitude  (p^, (R,T+) | 

Is  of  the  same  order  of  magnitude  as  the  populations  (p^^i  or  P3131)  in¬ 
duced  by  the  LICET  reaction.  Second,  states  2‘  and  3'  can  be  of  the  same 
or  opposite  parity  -  the  relative  parity  of  the  two  levels  determines 
which  part  of  the  interatomic  potential  (l.e.  dipole-dipole;  dipole- 
quadrupole,  etc.)  contributes  in  Eq.  (46).  Third,  th  LICET  cross  section 
is  significantly  enhanced  if  the  energy  of  the  intermediate  levels 

(e  or  e')  in  Eq.  (46)  is  such  as  to  lead  to  a  "small"  energy  denominator 
18 

in  Eq.  (49)  ;  in  that  case,  a  single  term  dominates  the  summation  in 

Eq.  (46).  The  '.-xlstence  of  nearly  resonant  intermediate  levels  has  nlaved 
a  key  role  in  all  experimental  observations  of  LICET  to  date. 

The  types  of  level  schemes  which  lead  to  near-resonant  enhancement 
have  been  discussed  in  a  previous  work.^  To  illustrate  this  feature,  we 
consider  the  level  schemes  depicted  in  Figs.  4  and  5.  In  Fig.  4,  contri¬ 
bution  to  Eq.  (46)  with  level  r’  as  intermediate  state  is  dominant.  One 
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can  view  the  LICET  process  as  a  collisional  exchange  which  takes  atom  A 
from  state  1  to  f  and  atom  A'  from  state  i'  to  r'  (the  fr'  state  Is 
virtual)  followed  by  the  field  E  taking  atom  A'  from  r'  to  2'  (leading 
to  the  amplitude  a^i)  —  field  E^  taking  atom  A'  from  r'  to  31 
(leading  to  the  amplitude  a^*)*  For  this  type  of  level  scheme,  levels 
2'  and  3*  must  be  of  the  same  parity.  The  value  of  produced  in 

this  LICET  reaction  may  be  calculated  from  Eqs.  (46 ) —(53)  and  Is  found 
to  be  of  order^ 


I  q^.^.T’X  -  K  ^  I  C  x'™ ' 

-  z  a  2. 

X  /  O  bA< 


(54) 


where 


*-2  A'  ^  tc-  ^ 

X'3-v 


(  ^  1  f.  ^  ^  ^  ^ 

<v  i  r;  \  a'  > 


(55) 


ur  is  the  most  probable  relative  speed,  b^f  Is  a  radius  which  characterizes 
resonant  broadening  of  the  r'-l’  transition,  and  b^  Is  a  radius  at  which 
the  collisional  level  shifting  operator  becomes  important  (as  such.  It 
Is  a  characteristic  radius  of  foreign  gas  broadening).  Values  of  b^,  In 
the  10X  to  2oX  range  are  typical  as  are  ratios  br,/bo  a  4. 
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It  would  appear  that  a  coherence  between  states  of  opposite 


parity  can  be  produced  by  the  level  scheme  shown  In  Fig.  5  In  which 
there  are  two  nearly  resonant  Intermediate  states  r'  and  r^  of  opposite 
parity.  Assume  that  states  1, 1' ,r^, 2 '  have  parity  "+"  and  states  2,r',3' 
have  parity  Then,  a  dipole-dipole  collislonal  interaction  followed 

by  field  E  acting  on  atom  A'  leads  to  a  final  state  amplitude  a^i  v^a 
the  pathway  21'  -*•  lr'  •+•  12'.  Similarly,  a  dipole-ijuadrupole  collislonal 

interaction  followed  by  field  acting  in  atom  A'  leads  to  a  final  state 
amplitude  a^«  via  the  pathway  21*  -*•  lr|  -*■  13'.  Since  the  collision 
operator  ~IL  contains  both  dipole-dipole  and  dipole-quadrupole  components, 
the  same  collision  creates  both  a^t  and  a^».  Consequently,  it  would 
appear  that  a  coherence  £>2' 3'  states  2’  and  3*  having  opposite 

parity,  can  be  produced  in  LICET.  The  order  of  magnitude  of  this  coherence 
is  given  by  an  equation  analogous  to  Eq.  (54). 

00  J 

Although  a  given  collision  produces  a  coherence  »  one 

finds  that  the  macroscopic  dipole  moment,  obtained  by  averaging  over 

all  possible  collision  orientations,  vanishes  if  states  2'  and  3' 

19  7 

are  of  opposite  parity.  In  some  sense,  one  can  view  the  production 
of  coherence  the  A'  atoms  as  a  four-wave  mixing  process; 

the  four  fields  are  the  two  laser  fields  and  the  dipole  and  quadrupole 
collislonal  interactions  acting  on  the  A'  atoms.  The  dipole  and 
quadrupole  collislonal  interactions  can  be  thought  of  a6  unpolarized 
"fields";^  for  isotropic  collisions,  these  fields  are  incident 
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from  all  directions  with  equal  probability  and  the  macroscopic  co- 
herence  vanishes.  However,  by  velocity  selecting  either  the 

A  or  A*  atoms  (e.g.  use  velocity  selective  excitation  of  the  A 
atoms,  use  a  beam  of  A  or  A*  atoms,  detect  only  those  A'  atoms  in 
a  given  velocity  subclass),  one  creates  an  anisotropic  distribution 
of  collision  orientations.  In  effect,  the  collisional  dipole  and 
quadrupole  "fields"  are  no  longer  incident  with  equal  probability 
from  all  directions  and  it  becomes  possible  to  create  a  nonvanishing 
macroscopic  A  detailed  calculation  of  p2'3»  be  given 

in  a  future  paper;  at  this  point,  however,  we  note  that  it  appears 
that  it  is  necessary  to  be  detuned  from  exact  resonance  for  one 
of  the  LICET  transitions  in  order  to  produce  a  nonvanishing  macro¬ 
scopic  coherence. 
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V.  Detection 


The  result  of  the  CARE  or  LICET  reaction  is  to  create  a  coherence 
in  one  of  the  atoms,  A  or  A'.  The  method  of  detection  depends  on  whether 
the  final  states  have  the  same  or  opposite  parity. 

A.  F-fnal  states  having  the  same  parity  following  the  CARE  or 
LICET  excitation. 

The  analysis  is  the  same  for  both  CARE  and  LICET;  to  be 
specific,  we  analyze  the  CARE  process.  Immediately  following  the  laser 
pulses  at  t  ■  T+  (which  is  now  arbitrarily  set  equal  to  zero) ,  a  coherence 


iVh-w.VR 


(5 


is  created.  Since  states  2  and  3  have  the  same  parity,  one  uses  an 

20 

interrogation  pulse  to  monitor  the  2-3  coherence.  Several  methods  are 
20 

available  ,  depending  on  the  relative  magnitudes  of  the  decay  rates  and 

the  k  and  k^  vectors.  We  shall  choose  one,  a  detection  scheme 

21 

similar  to  that  employed  in  the  tri-level  echo  ;  in  using  this  scheme, 
we  assume  that 

-  W,  \  u  >> 


where 


r. 


c 


is  the  total  rate  at  which  the  lj  coherence  decays. 


The  level  scheme  is  as  shown  in  Fig.  6.  At  time  after  the 


initial  laser  pulses,  a  third  laser  pulse 


L(k2'*  > 


in  near  resonance  with  the  2-4  transition,  is  incident  on  the  sample. 

The  pulse  duration  T  is  assumed  to  be  short  enough  so  that  all  relaxation 
processes  can  be  neglected  during  the  pulse  [i.e.  T  <<:  1; 

|k-k^|uT  <<  1;  k2uT  K<  1;  (f^  "  W42^  <<:  1*  etc.].  To  calculate  the 


response  of  the  system,  ve  set 


s,,  - 


(60a) 


in  anticipation  of  a  signal  being  generated  with  propagation  vector 

and  frequency  *2^.  In  addition,  we  write 

-c  c  kl  •*  -silt)  ( 

c  -  c’  .  e 

^  x  i  ^ 

t  [  (  w  -wl)  .*  -vp.+  -sO* 1  .  ( 

e' "  “ 


(60b) 


(60c) 


At  a  time  T~  just  before  the  interrogation  pulse  is  applied, 

*  0.  The  coherence  p 23  has  evolved  freely  to  a  value 
P23(R,0)  exp[  i(k-k1)  •  (R-vT^)  -Kfi-fi^Tj  -  3Ti ^ »  which*  together  with 
Eq.  (61b)f  implies  that  ^23^1^  *8  etlua*  t0 

- :  t  -  C  kx  -w.vwy.d 
■-  e  c 

i  [cif  -  (  T,  -.ya%T|  (( 
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Following  the  laser 


The  laser  field  couples  the  coherences  p2j  and  p^  • 
pulse,  acquires  a  value  proportional  to  p2j(R»T^);  with  the  neglect  of 
all  relaxation  during  the  laser  pulse  and  the  use  of  an  optimal  pulse 
(one  having  a  pulse  area  of  n/2),  the  value  of  p^  at  a  time  T*  immediately 
following  the  laser  pulse  is 

l  (62) 

For  t  >  T, ,  p  evolves  freely  according  to 
1  A3 

-  L  rA  ”  -^5)  T  L  W  1  £*3  -  (63) 

'  ^3 


Combining  Eqs.  (61)-(63)#  one  finds  that,  for  t  > 

„  -•»  L  kf  -  '  * 

--  Lqxtf,c)t 

l  ^  xi*.  _ ]  T,  (  n f  -  w-'-t  % )lt  * 1  >  ^ 


*  e 


*  e 


>  e. 


-  t  Cp 


(64) 


The  coherence  p. 

43 

phase-matching  condition 


generates  a  radiation  field  provided  that  the 

-9 

-  K,-k,'k 


(65) 
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along  with  the  subsidiary  condition 


< 


C-H-t  '  )  L  /c 


*  1 


(66) 


Is  obeyed  (L  *  sample  length).  These  phase-matching  conditions  can  be 
achieved  by  taking  copropagating  waves  with  (flj-fl)=(i)  ^  end  ^2=W42*  ABSumin8 
the  phase-matching  conditions  to  hold  and  that  the  sample  is  optically  thin, 
it  is  easy  to  show  that  the  power  density  exiting  the  sample  is  given  by 


v.o 


s  - 


hi 


|<  ^  I  f*  I  I  ,7C  )  L  1  ^ 


(fc  v’  t)> 


(6?) 


where  N  is  the  atom  A  density  and  the  average  is  over  the  atom  A  velocity 

A 

distribution  and  any  changes  in  caused  by  the  lack  of  spatial  co¬ 
herence  of  the  laser  fields.  On  integrating  over  velocities  for  times 

(t-T.)  >  l/k,u,  one  finds  a  negligibly  small  p  except  when 
it  43 

k>T.)  Hk-iOT,  '  0  •  (68) 


For  copropagating  waves  satisfying  Lq.  (65),  this  condition  implies  that  an 

echo  can  be  produced  at  a  time  t  given  by 

€1 


t 


(69) 


An  echo  is  produced  provided  that  k^  >  k,  an  inequality  that  holds  for  the 
system  we  have  chosen.  Combining  Eqs.  (64)— (68) ,  we  find  that  the  maximum 
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value  for  the  echo  power  density  exiting  the  sample  is 


K  /  *  ,  i 

C-  (t  -  t  _  l  <  I  P»  i  *  >  t  (irrc  )  Wf  L 

-  *\iz\  ts  +  r^et^rj] 

«.  e 


\  <  ?„l*.  <=)>'" 


(70) 


Asbuoins  that  j  < p „ ^ (R , 0 ) > |  can  reach  values  of  order  0.01,  one  finds  power 

2  15  3 

densities  of  order  l.OW/cm  for  active  atom  densities  N.  =  10  atoms/cm  . 

A 


b.  Final  states  having  opposite  parity. 

For  the  LICET  reaction,  it  is  possible  to  create  an  optical 
coherence  p ^ »  in  which  states  2’  and  3'  have  opposite  parity.  This  co¬ 
herence  can  be  monitored  by  detecting  the  f ree-induction  decay  signal 
emitted  by  the  sample.  At  time  T+  ■  0  immediately  following  the  LICET  ex¬ 
citation,  a  coherence 


-> 
R  . 


-  + 

i 


--  o)  e. 


x(>'  **  tV 


ft 


(71) 


i6  created.  This  coherence  then  evolves  freely  so  that  at  a  time  t  >  0,  one 
finds  phase  matching  can  be  achieved  with  copropagating  fields  having 
k^-k  =  ^ «  *  Under  phase-matching  conditions,  the  power  density  exiting 

the  sample  is  given  by  an  equation  analogous  to  (67).  On  performing  the 
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necessary  velocity  integration  assuming  (k^-k)u'  >>  • 3 *  (u'  “  atom  A' 

most  probable  speed),  one  may  obtain 


S  t  ^ 


|  <  »'  1  fK'  1  «-'  >  l‘  (‘>‘■1  u.- k  )’L‘ 

|  <r  42,,-  <~Kc  )>  I 


7~ 


(72) 


where  N  ,  is  the  A'-atom  density.  Although  the  value  for  the  coherence 

A 

I ^2  *  3 ’ (R»  0) i  is  smaller  than  that  in  the  case  when  levels  2'  and  3'  have 
the  same  parity,  (see  discussion  of  Sec.  IV),  it  would  appear  that  there 
is  sufficient  signal  strength  to  detect  the  LICET  induced  coherent  emission. 
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VI.  Conclusion 


A  method  for  producing  electronic  state  coherences  using  Combined 
Field-Collisional  Transitions (CFCT)  has  been  outlined.  Such  coherences 
can  be  produced  using  either  Collisionally-Aided  Radiative  Excitation 
(CARE)  or  Laser-Induced  Collisional  Excitation  Transfer  (LICET).  The 
final  state  coherence  may  be  created  between  levels  having  the  same  or 
opposite  parity.  For  final  states  of  the  same  parity,  the  coherence  can 
be  "stored";  at  some  later  time,  an  interrogation  pulse  can  be  used  to 
trigger  the  emission  of  dipole  radiation.  For  final  states  of  opposite 
parity,  the  system  can  radiate  immediately  following  the  LICET  excitation. 
In  this  way,  one  can  use  LICET  for  either  sum  or  difference  frequency 
generation. 

As  noted  above,  some  additional  feature  must  be  added  to  LICET 
(i.e.  velocity  selection  of  either  the  A  or  A'  atoms)  to  produce  a  final 
state  coherences  between  states  of  opposite  parity.  Another  way  tb  achieve 
a  final  state  coherence  between  .  tates  of  opposite  parity  using  either 
CARE  or  LICET  is  illustrated  in  Fig.  7.  A  coherence  is  created  between 
two  states  of  atom  A  (i  and  i^)  by  a  CW  or  pulsed  field.  From  this  Initial 
state,  a  CARE  reaction  produces  a  coherence  p^2  (fig*  7a)  while  a  LICET 
reaction  produces  a  coherence  p2'$'  7b)  between  states  of  opposite 

parity.  This  technique  could  be  used  to  produce  radiation  at  the  sum 
frequency  of  the  three  laser  fields. 

In  order  to  produce  coherence  using  CFCT,  the  laser  fields  must  be 
temporally  and  spatially  coherent.  The  relative  phase  of  the  two  laser 
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fields  cannot  vary  significantly  in  times  less  than  the  pulse  duration  or 
in  distances  less  than  a  wavelength  for  the  coherence  to  be  produced. 

Since  the  production  of  tri-level  echoes  requires  the  same  coherence 
properties  and  since  tri-level  echoes  are  readily  observed  with  non- 
mode-locked  lasers  ,  it  would  appear  that  the  coherence  criteria  for 
creating  CFCT  induced  electronic  state  coherences  can  be  achieved. 
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Figure  Captions 


1.  CARE  reaction  +  A|,  +  ftft  A(12)  +  A^ ,  [the  notation  A(1J)  Is  used 
to  Indicate  a  coherence  between  states  1  and  j].  The  A  and  A'  atoms 
collide  in  the  presence  of  a  pulsed  laser  field  whose  frequency  is 
represented  by  an  arrow  in  the  figure.  As  noted  in  the  text,  the 
coherence  p  vanishes  on  averaging  over  times  during  the  laser  pulse 
at  which  the  collision  may  occur. 

2.  CARE  reaction  A^  +  A£ ,  +  Tift  ♦  fift^  -*■  A(23)  +  A|,.  A  coherence  is 
created  if  ft^  -  ft  - 

3.  LICET  reaction  A£  +  Aj ,  +  Tift  +  -*•  Ax  +  A'  (2 ' 3' )  .  A  coherence  p2,3, 

is  created  if  -  ft  ■  3a)  or  ft^  +  ft  =  3b). 

4.  LICET  reaction  A^  +  A^(  +  fift  +  -fift^  A^  +  A,(2'3')  with  a  nearly 
resonant  intermediate  state.  The  colllslonal  interaction  creates  a 
virtual  state  with  energy  close  to  that  of  level  r'  and  the  field 
interactions  complete  the  LICET  process.  States  2'  and  3'  must  have 
the  same  parity, 

5.  A  LICET  reaction  similar  to  that  shown  in  Fig.  4,  but  one  for  which  there  are 

two  nearly  resonant  Intermediate  states  of  opposite  parity  and  for  which 
states  2'  and  3'  are  of  opposite  parity.  The  r'  and  r^  levels  en¬ 
hance  the  excitation  probability  for  states  2'  and  3',  respectively. 
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Figure  Captions  -  Con't. 


6.  A  scheme  for  detecting  the  coherence  created  by  a  CARE  reaction 

at  t  »  0.  At  time  a  laser  pulse  of  frequency  fij  =  *8  aPP^^e^ 

to  the  atoms.  At  some  time  later,  an  echo  signal  can  be  generated 
at  frequency 

7.  A  method  for  producing  electronic  state  coherences  between  states  of 
opposite  parity  by  starting  with  atom  A  prepared  in  a  coherent 
superposition  of  states  1  and  i^. 

(7a)  CARE  reaction  Ad^)  +  A^ ,  +  -*•  A(23)  +  A^, ; 

(Q  +  fil  +  =  W32  ‘ 

(7b)  LICET  reaction  A(ii1)  +  Aj ,  +  ^  +  A'  (2 ' 3' ) ; 

+  “l  +  Uii,i)  =  W3 ' 2 '  * 
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The  problem  of  calculating  transition  probabilities  in  two-level  systems  is  studied  in  the 
limit  where  the  detuning  is  large  compared  to  the  inverse  duration  of  the  interaction.  Cou¬ 
pling  potentials  whose  Fourier  transforms  * '( I  are  of  the  form  / ire  le  for  large  fre¬ 

quencies  give  rise  to  solutions  which  may  be  classified  into  families  according  to  the  form  of 
Within  each  family  transition  probabilities  may  be  calculated  from  formulas  that 
differ  only  in  the  numerical  value  of  a  scaling  parameter.  In  cases  where  the  coupling  func¬ 
tion  has  a  pole  in  the  complex  time  plane,  the  families  are  identified  with  the  order  of  this 
singularity.  In  particular,  for  poles  of  first  order,  a  connection  with  the  Rosen-Zener  solu¬ 
tion  can  be  made.  The  analysts  is  performed  via  high-order  perturbation  expansions  which 
are  shown  to  always  converge  for  two-level  systems  driven  by  coupling  potentials  of  finite 
pulse  area. 


I  INTRODUCTION 

In  many  areas  of  physics,  one  encounters  prob¬ 
lems  involving  two  states  of  a  quantum-mechanical 
system  coupled  by  a  time-dependent  potential.1-11’ 
In  the  interaction  representation,  the  equations  of 
motion  for  a,  and  a:,  the  probability  amplitudes  of 
levels  1  and  2,  are  of  the  form 

to i  =  FiMe'“"cM  ,  (la) 

id :  =  Ft  Me  ~,u,a ,  ,  (lb) 

where  w  is  the  frequency  separation  of  the  states  and 
FI M  is  the  coupling  potential.  Decay  effects  are 
neglected  in  Eqs.  (1)  (and  throughout  this  paper), 
and  we  work  in  a  system  of  units  in  which  /}  =  1. 

Equations  of  this  type  arise  in  many  semiclassical 
problems.  A  problem  of  current  interest  to  which 
they  apply  is  the  coupling  of  two  levels  of  an  atom 
by  a  laser  pulse  that  has  a  temporal  width  which  is 
small  compared  to  the  natural  lifetimes  of  the  levels. 
The  pulse  F(U  is  of  the  form 

F( t )  =  2.4 ( i icoslh  ,  (2) 

where  11  is  the  central  frequency  of  the  pulse,  and 
2,4  it)  is  the  envelope  function  of  its  amplitude.  As¬ 
suming  that  il  /til  4  (i>  )<-"!,  one  can  recast 
Eqs.  (1)  in  terms  of  A,  the  detuning  of  the  pulse 
from  resonance  (rotating-wave  approximation),  as 

id  |  =5.-1  1  /  h,lAla :  ,  (3a) 

ta ,  -  A  i  Me  lA'a|  .  1 3b) 

Equations  (3)  or  i|)  are  deceptively  simple  in 
form,  and  one  might,  at  first  glance,  believe  that  the 
system  must  he  completely  understood,  so  that  noth- 
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ing  remains  to  be  investigated  about  the  equations  or 
their  solutions.  Actually,  there  is  very  little  known 
about  the  overall  qualitative  nature  of  the  solutions 
to  Eqs.  (3)  for  arbitrary  Ail).  Apart  from  any  in¬ 
trinsic  interest  one  might  have  in  the  dynamics  of 
two-level  systems,  such  information  could  be  useful, 
for  example,  in  applications  where  one  wished  to 
choose  the  pulse  shape  to  maximize  the  excitation 
probability  for  a  given  detuning  A. 

To  appreciate  that  our  assertion  concerning  the 
lack  of  knowledge  about  the  behavior  of  systems 
described  by  Eqs.  (3)  is  valid,  one  need  only  recog¬ 
nize  that  the  answer  to  the  following  question  is  not 
known  in  general:  Starting  with  initial  conditions 
a,(  —  x  )  =  1  and  u:(  —  x  >  =  0,  how  does  the  proba¬ 
bility  amplitude  a2it)  depend  qualitatively  on  the 
pulse  area  S,  defined  by 

S  =  f  *  A  ( t'Jt  , 

on  the  detuning,  and  on  the  shape  of  the  envelope 
function  Air)?  A  response  to  this  query  can  be 
made  for  a  limited  number  of  cases.  Analytic  solu¬ 
tions  are  available  if  A  it)  belongs  to  a  class  of  func¬ 
tions'  (including  the  hyperbolic  secant  of  Rosen  and 
Zener-1)  mappable  into  the  hvpergeometric  equa¬ 
tion,  or  if g- 10 

A  U)  =  ( const  )exp(  -a  ]  f  I  )  , 

or  if  A  it)  is  a  step  function  (Rabi  problem),  or  if  the 
detuning  is  zero.  (Kaplan  has  also  considered  cases 
where  the  detuning  varies  as  prescribed  functions  of 
the  amplitude  and  obtained  closed-form  expres¬ 
sions.)  In  addition,  there  are  approximate  solutions 
available  in  adiabatic4  or  perturbative  limits.  Yet, 
there  remains  a  wide  range  of  parameters  and  pulse 
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■Tapes  I\ii  wbuh  an  answer  lo  the  basic  quc-tion 
e.mili'l  he  prov  nled 

In  this  papci .  wi1  -^1  < .t  1 1  exan  ,nv  .hi-  solutions  to 
Fqs  3  in  the  limit  where  the  produ  t  of  the  detun¬ 
ing  A  .m J  the  characteristic  pulse  ihiralii'U  r  has 
a  mucnitude  greatly  m  excess  of  units.  In  other 
w  ids.  »e  are  assuming  that  the  pulse  does  not  pos¬ 
sess  the  appropriate  Fourier  components  to  sigmfi- 
v.mtlv  compensate  for  the  detuning.  Ill  conse¬ 
quence.  the  transition  probability  uo  x  1  '  will  ai¬ 
rs  a\s  be  set'  small  hut  still  great  enough  to  he  ex¬ 
perimentally  measurable  in  atomic  sapors  of  densi¬ 
ties  -  atoms,  cm  '1.  We  note  that  numerical 
solutions  of  Fqs.  oh  in  this  detuning  range  may  be 
possible  hut  are  sets  costly  in  computer  time  and 
plagued  with  technical  difficulties. 

For  the  case  At  .  .  1,  we  shall  establish  the  fol¬ 
lowing  results  1 1 1  I.oss-order  perturbative  approxi¬ 
mations  lor  ,/ •  y  are  not  valid  for  arbitrate  pulse 
area  A.  despite  the  fact  that  a t '  1  for  all 

time  2  Ail  iterative  solution  to  Fqs.  i|i  always 
conserges  (or  well-behaved  envelope  functions.  (3) 
Asymptotic  solutions  for  u  :l ft,  I  finite,  may  be  easi¬ 
ly  found,  but  expressions  for  up  x  >  are  difficult  to 
obtain.  41  Asymptotic  solutions  for  u : I  -x  i  can  be 
obtained  for  a  limited  class  of  pulse-envelope  func¬ 
tions  using  contour  integration  techniques.  This  is  a 
broader  set  than  that  for  which  exact  solutions  are 
known  51  The  asymptotic  dependence  of  a,l  x  i 
depends  critically  on  the  nature  of  the  singularities 
of  the  pulse  -  envelope  function  .4  in,  analytically 
continued  into  the  complex  plane,  t ft i  If  two  pulse 
functions  have  the  same  Fourier  transforms  in  the 
limn  of  large  frequencies  and  if  the  dominant 
dependence  of  the  transform  is  ;m  exponential  decay 
m  the  frequency,  then  the  asymptotic  forms  of  the 
solutions  a p  -x  1  for  these  functions  in  the  limit  of 
large  A  are  simply  related.  In  this  paper,  we  ad¬ 
dress  points  1 1 1,  (2),  (3i,  and  (6);  methods  for  actually 
obtaining  asymptotic  solutions  [points  (4i  and  (5t] 
will  be  discussed  in  a  future  article.  In  the  present 
discussion,  the  initial  conditions  are  taken  as 
u  |  x  i  I  and  a 2 1  -  -x  )  -  0. 

II  ASYMPTOTIC  SOLUTIONS 

As  we  have  indicated,  the  Rosen-Zener  (hyper¬ 
bolic  secant  coupling  pulse)  problem  is  one  of  the 
few  for  which  exact  solutions  are  known.  In  this 
case,  a  simple  expression  gives  the  transition  ampli¬ 
tude  as  a  function  of  detuning  ;md  area  for  all  values 
of  these  para,  .cters.  Naturally,  since  tins  formula 


w  here  t  is  the  Fourier  t  ransf  i  mn  of  1  1  l  ,  is  exact ,  it 
is  valid  in  the  special  case  ol  the  asvmptoli.  limit 


We  shall  show  that  there  is  an  culiic  class  o| 
pulses  for  which  the  asymptotic  tiaiisition  umpli 
Hide,  as  a  f  unction  of  S  and  A.  in.iv  tie  w  lit  ten  by  in 
sped  ion  once  the  Rosen  -Zenci  pioblem  has  been 
solved.  We  shall  also  demonstrate  that  ibeic  are 
other  classes  of  pulses  whose  solutions  as  ;  .  ,  ate 

unrelated  to  Rosen-Zener  but  ale  coimeeled  to  each 
other  in  tfie  sense  that  once  one  has  been  solved,  the 
solutions  for  tile  enure  class  may  be  obtained  by  in¬ 
spection. 

The  existence  of  these  related  solutions  will  he  es¬ 
tablished  via  term-by-term  comparison  of  nth-order 
perturbation  expansions  which,  under  very  general 
conditions,  arc  convergent  in  two-level  problems  tsec 
the  Appendix).  With  suitable  sealing  of  the  cou¬ 
pling  strengths,  the  series  for  different  members  of 
particular  classes  will  be  seen  to  be  identical  in  the 
limit  of  large  detunings. 

The  particular  potentials  analyzed  in  this  paper 
are  .-1 '  / 1  whose  Fourier  transforms  for  large  <■>  as¬ 
sume  the  form  ploilexp i  —  boj  i,  where/)  is  slowly 
varying  in  a  frequency  interval  b  \  and  b  is  a 
constant.  It  is  convenient  to  make  a  variable 
change,  such  that  v  =  b  uj  and  .v  -  t  /  b  .  Conse¬ 
quently,  the  exponential  decay  factor  in  the  Fourier 
transform  becomes  expi  v  i  and  the  equations  of 
motion  transform  to 


ia  i  -  flf1  x  k’lr,xa ;  ,  <3a' i 

iii  2  -  ft f 1  v  'e  ",xu  i  ,  1 3b'  i 


where  a  />A  and  where  the  dot  now  signifies 
differentiation  with  respect  to  .v.  The  quantity  f). 
previously  designated  as  .V.  is  the  pulse  area.  The  re¬ 
duced  potential  function  / i.vl  is  defined  such  that 

j~  ftxkix  -  1  . 


The  pulse  area  is  invariant  under  the  indicated 
change  of  variable.  One  may  also  w  rite  Fqs.  <  3 )  as  a 
pair  of  uncoupled  second-order  equations 


4  i  a 


d,  0  , 

ift  •  ///«■  0  . 


oa> 


iS|n 


There  arc  two  aspects  to  the  solutions  ol  Fqs.  3 
or  (5).  T  hese  tire  the  ealeuhiiions  of  the  amplitudes 
at  finite  and  infinite  times,  respeetr  ely.  I  he  former 
are  of  interest  il  the  transient  solutions  are  to  be 
useil  as  inputs  to  other  proK  s.  such  as  multipho 
ton  ionization.1 !  while  the  lat'  with  which  we  are 
mainly  concerned  here,  gives  the  transition  ampli¬ 
tude  u  /  '.  I  he  two  temporal  regimes  diflci  great 
ly  m  the  methods  that  must  be  used  tn  perform  ac 
curate  e.ih  illations 
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One  may  write  the  solutions  to  Eqs.  O'  as  pertur¬ 
bation  series  in  the  usual  fashion,  noting  that  only 
even  orders  enter  the  expression  for  a , .  while  only 
odd  orders  appear  in  the  formuhi  for  </;.  The  ex¬ 
pansion  for  ti;i  •  x  1  is 

a:  -  i  2  a‘:lk''-p2k  +  ' i 

A  -i) 

where 

■  j*  f  i  rr  . 

a :  =  J  J 1  -v  |  )e  (/.k  j 

■  II  J  x  J'x,w  -dxj  . 


In  the  Appendix,  it  is  shown  that  this  series  eon- 
verges  for  all  finite  pulse  areas 

Tor  the  remainder  of  the  paper  we  will  restriet 
ourselves  to  the  ease  of  pulses  that  are  symmetric  in 
time  and  where  a  -si,  the  adiabatic  or  asymp¬ 
totic  limit.  The  Fourier  transform  will  be  sym¬ 
metric  in  v.  We  shall  begin  by  comparing  the  finite 
and  infinite  time  solutions  of  the  Rosen-Zener  prob¬ 
lem,  which  exemplify  relevant  properties  of  transi¬ 
tion  amplitudes  induced  by  smooth  pulses. 

With  initial  conditions  a  |(— x)-  I  and 
till  -  x)  =  0  with  a  pulse-envelope  function 
/(.t)  =  sechi  ~.x/2)/2,  Rosen  and  Zener1'  obtained 
an  analytic  solution  to  Eqs.  (3'f  of  the  form 


1 6a  I 
I6b! 

(6b'l 


and  /  designates  the  hy  pergeometnc  function.  The  form  of  a:  given  by  Eq.  (6b)  is  valid  for  all  x,  while  that 
given  hv  I  q  6b  holds  only  for  finite  v,  unless  []  corresponds  to  an  eigenvalue,  a  pulse  area  for  which 
j  ■  /  vanishes  '  We  recall  that  u- 1  x  i,  the  transition  amplitude  for  the  Rosen-Zener  problem,  is  given  by 
Eq  4 

We  m.iv  obtain  the  finite  time  solution  by  explicitly  expanding  the  2F ,  function  of  Eq.  (6b') 


For  large  <r,  it  is  sufficient  to  retain  the  leading  term 


This  is  equivalent  to  first-order  perturbation  theory 
in  the  adiabatic  limit 

u\'  if)  f'  fix'  ie  lax'dx ' 

o/'-V  ,  mx 
- 1>  «' 
a 

where  subsequent  parts  integrations  arc  neglected, 
since  they  are  O'  l /it"),  n  •  1.  We  immediately  see 
that  this  sequence  of  parts  integrations  is  unsuitable 
for  calculating  a- 1  /  ',  since  each  term  separatelv 


vanishes  when  ,x  — •  x .  Even  including  the  third- 
and  higher-order  terms  in  the  perturbation  series  via 
analogous  sequences  of  parts  integrations  does  not 
enable  one  to  obtain  a  nonzero  amplitude  as  i  *  x  . 
Consequently,  other  methods  are  necessary  to  calcu¬ 
late  (l.i  oc  I. 

It  is  clear  from  the  preceding  paragraph  that  for 
large  enough  or,  first-order  perturbation  theory  is  a 
sufficiently  accurate  approximation  for  most  pur¬ 
poses,  provided  x  is  finite.  For  infinite  times,  not 
only  does  the  adiabatic  sequence  of  parts  integra¬ 
tions  lead  to  an  incorrect  «-(  x  1.  but  even  an  exact 
evaluation  of  the  first-order  integral  may  he  insuffi¬ 
cient.  Ellis  is  typified  by  the  exact  Rosen-Zener  am¬ 
plitude,  Eq.  i4h  in  which  the  factor  sin fi  does  not 
reduce  to  its  first-order  limit  of  fi  unless  fi  is 
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small  compared  to  unity.  This  failure  of  the  first- 
order  theory  occurs  no  matter  how  large  the  detun¬ 
ing  becomes.  One  must  retain  enough  terms  in  the 
perturbation  expansion  to  accurately  represent  the 
sine  function.  Thus  for  t fie  Rosen-Zener  pulse,  if 
the  coupling  is  great  enough  so  that  saturation  ef¬ 
fects  would  appear  at  resonance,  simple  first-order 
theories  cannot  be  used  for  a  nonresonant  pulse  of 
the  same  strength.  As  we  shall  see,  other  smooth 
pulses  also  possess  this  "saturation  memory."  In 
fact,  in  some  cases,  a  higher-order  theory  is  neces¬ 
sary  off  resonance  even  for  a  case  where  a  first-order 
theory  would  suffice  at  resonance.  This  is  exempli¬ 
fied  by  the  formulas  of  Eqs.  (9)  below. 

Since  each  coupling  function  fix)  is  different,  one 
might  be  led  to  believe  that  separate  calculations 
must  be  performed  for  each  individual  case.  For¬ 
tunately,  as  we  have  stated  earlier,  there  prove  to  be 
classes  of  pulses  where,  if  one  knows  the  functional 
dependence  of  the  asymptotic  transition  amplitude 
on  «  and  ft  for  one  member  of  the  class,  one  knows 
it  for  all  members  of  the  class,  although  the  actual 
time  dependence  of  the  potentials  may  be  drastically 
different.  What  is  significant  is  that  their  Fourier 
transforms  assume  the  same  form  as  a  — •  x  . 

When  Rosen  and  Zener  deduced  Eq.  (4),  they  sug¬ 
gested  that  similar  formulas  might  hold  for  other 


smooth  pulses  ~  This  conjecture  proves  not  to  hold 
in  general.  It  is  manifestly  false  for  asymmetric 
pulses  and  is  not  even  valid  for  all  symmetric 
pulses.'-''  What  we  shall  show  is  that  a  kind  of 
Rosen-Zener  conjecture  does  apply  at  large  detun¬ 
ings  for  pulses  in  which  fix'  has  simple  poles  at 
x  i.  This  law  does  not  apply  to  pulses  which  have 
higher-order  poles  at  tilts  point,  although  scaling 
laws  for  these  do  exist,  different  for  each  order. 

The  following  theorem  will  be  established.  Let 
two  coupling  pulses  fix)  and  have 

Fourier  transforms  / tv)  and  /,,(vi.  The  Fourier 
transforms  of  both  approach,  for  large  values  of  the 
argument,  the  same  asymptotic  form  faiv  i.  1  f  fa  is 
of  the  form  Si  v  k>  ~  v  ,  where  Si  v  )  is  a  slowly  vary¬ 
ing  function  of  v,  then  the  asymptotic  transition  am¬ 
plitudes  generated  by  the  two  pulses  will  be  the 
same,  provided  that  the  pulse  areas  are  both  finite. 
A  sufficient  condition  for  the  indicated  asymptotic- 
behavior  of  the  Fourier  transforms  is  that  they  be 
equal,  for  large  v,  to  a  contour  integration  whose 
value  is  given  by  the  product  of  the  residue  at  .v  / 
and  the  usual  Cauchy  factor  2-i.  If  two  such  pulses 
are  to  have  the  same  <A(vi,  they  must  possess  poles 
of  the  same  order  at  x  —i. 

The  contribution  of  order  ilk  -  1 1  to  the  transi¬ 
tion  amplitude  may  be  rewritten  slightly  . 


**  —  nc 


2k  +  1 


*,)<?  'aX'dx]  ri  lim  f 

j  A;  -0  J 


c,  ,  It'  -  I 
Jix,  )e 


-l-b 


dx, 


The  factors  e  ’  1  do  not  affect  the  integrals.  They  are  used  to  remove  ambiguities  as  x  •  v.  m  the  treat¬ 
ment  below,  where  we  express  the  amplitude  in  terms  of  integrals  in  the  frequency  domain.  The  limits  /  *0 

are  to  be  taken  before  the  xt  integration  is  performed.  Expressing  each  <x;>,  j  >  2.  in  terms  of  its  Fourier 
transform,  we  find 


■Ik  +  I 


<2jt) 


1  r  * 

L  f  /(*  1 

rr  *  u  x 


2*  +  l 


dx'  n  ,iim„  //,  dxi  f  \,f{yj k‘l 


1  Ja  - 1 A 


I  tv, 


By  working  in  the  frequency  domain,  we  shall  be  able  to  examine  the  structure  of  the  integrals  for  a ;k  * 1  and 
establish  that  the  contribution  from  regions  where  the  asymptotic  form  of  /  is  not  valid  is  lower  by  Oi  1/ui 
than  the  contributions  from  regions  where  it  is  valid. 

The  integrals  over  the  x;  are  trivial  to  perform.  We  obtain 


'  2k  +  I 


lim 

i  -o  ( 2xr  >‘ 


r  00 

2k  +  1 

/  d  V2 

*  —  00 

'  dv2k  +  if 

1  vj  ~ a 

j=2 

2k  +  t 

I 


/  =  :*  +  n-j 


fi  V;  > 
[v,  +  t  - 


1 )'« 


We  now  proceed  to  determine  the  asymptotic  form  of  these  amplitudes.  The  analysis  is  easiest  to  follow  for 
the  third-order  contribution  ay1',  but  exactly  the  same  reasoning  and  conclusions  will  apply  for  the  higher- 
order  terms.  (The  theorem  is  true  by  inspection  in  first  order,  since  that  contribution  is,  apart  from  a  constant 
multiplier,  just  the  Fourier  transform  itself.  Thus  if  two  coupling  functions  have  Fourier  transforms  of  the 
same  asymptotic  form,  their  first-order  transition  amplitudes  scale  the  same  way  with  ft  and  a.)  The  leading 
nontrivial  term  is  o  V  Changing  the  dummy  variable  V|  to  v,,  we  find 


lim 

k  .a  i  2~ 


j  i  v,  1/ 1  v;  I/1  v,  *  v ,  -  a  )dv|  dv- 
t  »-|  n  //.  'I  v,  *  i'|  iA  i 
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where,  without  loss  of  generality,  all  A;  and  sums  of  A;  have  been  replaced  by  the  single  infinitesimal  a.  It  is 
convenient  to  make  the  change  of  variable  v,  ypr.  One  finds 


i  rx  rx  /l<ri’t 

hm  -  1 

1  f'ti.V;  './(otyi 

•  yq  -  1  i  h/l1 1  dy  - 

K  -n  V  2,7  J  J  * 

yq  1  /Aliy 

*  V:  -  if- 1 

i  p  r*  rx  S'“-v>V 

i  <nq  '/(«iy  |  *■ 

tq  -  1 1  )dy  |  dy : 

v  2  7?  ^  X  ~  x 

ty j  1  a y ;  s 

yq  i 

. 

r  1 

/  '  a  i|  /  wry  ,  i]‘ 

[/'  oyq  i  ]y/  i  a  1  1 

•  t,r  Inn  <1  v 

*- 

*  .  ^ 

1  -  l  /  A 

-y:  1  - / a  j 

where  /’  indicates  that  the  mtc  rand  excludes  infinitesimal  regions  near y,  —  -  yq  andy’i  =  1.  We  may  formally 
integrate  the  last  two  terms  II  !  is  factored  from  the  second  of  the  two  integrals,  they  combine  to  become 


r  y  1 

1,7  lull  I  (Jv  ■  /  '  It  l[  /  ‘  <<  V  • 

*.  ■■■'  J  *  '  "  '  '  i  v-  I  A 


It  is  immediately  obvious  that  if  these  are  parti¬ 
tioned  according  to  the  rule 


lim  f  6,x'd*  />  f  t,T, 

i  •  ■  X  V  '(  J  .X  X 


<3  1  x  'd.x 


<3  -x  „  ,  , 


the  principal  value  contributions  exactly  cancel, 
while  the  17  terms  are  proportional  to  e  and  are 
exponentially  small  compared  to  a  ,1  ,  which  decays 
only  like  e  feints  proportional  to  exponentials 
which  decay  more  rapidly  than  e  11  do  not  contri¬ 
bute  to  the  asy  mptotic  form. 

VV'e  now  proceed  to  examine  the  remaining  contri¬ 
butions  '  i  a  V  ,  where  it  is  again  understood  that  the 
small  re,  he  neighborhood  of  y ,  -  -  y,  and 

.V’l  =  1  are  i  from  the  integrals.  For  all  re¬ 
gions  i  are  y  <  a /a  .  where  a 

is  a  numbv  cr  unity,  fitly)  may  be  replaced 

by  its  asymptotic  torn)  fa(ny >.  Thus  for  the  entire 
V’  i  —  plane,  except  where  y  \  -  0,  yq-()  'but  not 
both  simultaneously)  and  yq  yq  "  1,  the  numerator 
of  the  integrand  is  well  represented  by  its  asymptotic 
form.  Furthermore,  since  at  most  one  of  the  three 
Fourier-transform  factors  departs  from  its  asymp¬ 
totic  form  in  any  given  region  of  space,  the  area  in 
the  yq  yq  plane  over  which  one  of  the  /  both 
departs  from  its  asymptotic  form  and  decays  no 
more  rapidly  than  e  a  is  ()i]/ct).  It  is.  of  course, 
implicitly  assumed  that  the  exact  and  asymptotic 
forms  of  the  Fourier  transforms  remain  bounded  as 
their  arguments  approach  zero.  For  the  former,  this 
is  equivalent  to  the  requirement,  which  we  have  al¬ 
ready  stated,  that  fi  be  finite. 

Now  consider  that  portion  of  the  yq  v;  plane 
where  all  factors  in  the  numerator  are  well  approxi¬ 
mated  by  their  asymptotic  forms.  Examine  m  par¬ 
ticular  the  exponential  decay  factors 


The  only  portion  of  the  plane  where  the  combined 
effect  of  the  exponential  factors  leads  to  an  overall 
decay  that  is  not  faster  than  e~a  is  the  range 
()<>'i<  1,  Oxy,  <  I  -yq.  The  integrand  does  not 
change  sign  in  this  portion  of  yq  —  yq  space,  which 
encompasses  an  area  -  t  (to  be  compared  with  the 
area  of  order  1  /«  which  was  found  for  the 
nonasymptotic  contribution).  Note  that  there  is  no 
portion  of  the  plane  in  which  the  integrand  decays 
more  slowly  than  e~a.  Thus  the  nonasymptotic  in¬ 
tegrand  contribution  is  O  (l/o)  compared  to  that  of 
the  asymptotic  integrand.  Similar  considerations 
enable  one  to  deduce  that  one  may  also  replace  the 
Fourier  transforms  in  the  higher-order  integrals  by 
their  asymptotic  forms.  We  thus  conclude  that  if 
the  time  dependences  of  two  coupling  functions  are 
such  that  the  asymptotic  forms  of  their  Fourier 
transforms  are  identical  and  of  the  indicated  form, 
the  large-detuning  transition  amplitudes  are  the 
same. 

As  we  have  indicated,  a  sufficient  condition  that 
two  pulses  have  the  same  a- t  x  i  for  large  a  is  that 
both  asymptotic  Fourier  transforms  be  equal  to  con¬ 
tour  integrations  given  by  (27i)[Res(x  =/')].  We 
compare  the  hypeibolic  secant  of  Rosen  and 
Zener,  /  =  ysecht 7.x  /2),  with  the  Lorentzian 
/  -  1 1/7 )( 1  v.x:r  '.  The  corresponding  .4  t.x) 
—  fif  t.x)  are 

th  ,  , 

A,  t.x)  1  I  +  ATI  '  . 


,  Pit  ,  7.X 

■hi  "  i  ,  sech  ^  . 

The  transforms  for  both  may  be  calculated  via  con¬ 
tour  integrations  The  I  orent/ian  case  is  trivial  and 
applies  to  all  v,  not  gist  large  frequencies.  We 
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choose  a  contour  that  runs  along  the  real  axis  from 
—  R  to  +  R  and  is  closed  by  a  semicircle  in  the 
upper  half  plane.  The  contribution  to  the  contour 
integral  from  the  arc  vanishes  as  R  *  x ,  so  that  the 
Fourier  transform  is  identical  to  the  contour  in¬ 
tegral,  whose  value  is  determined  by  the  residue  at 
the  simple  pole  at  x  =/.  The  result  is 

fti 

A,  =  e  '  .  (7a) 

L  V  2  77 

For  the  hyperbolic  secant  we  choose  a  rectangular 
contour  which  runs  from  ~R  to  +R  along  the  real 
axis,  and  is  continued  by  rectangular  segments 
parallel  to  the  imaginary  axis  from  the  points 
(±R, 0)  to  the  points  (±R,2i),  and  is  closed  by  a  line 
parallel  to  the  real  axis  which  runs  from  (R, 2i)  to 
(  —  R,2i).  The  two  vertical  segments  give  vanishing 
contributions  as  R— *&>,  and  the  horizontal  segment 
off  the  real  axis  goes  exponentially  to  zero  compared 
to  the  segment  along  the  real  axis  as  v— «•  oc .  Thus 
for  the  hyperbolic  secant,  the  Fourier  transform  is 
identical  to  that  of  the  Lorentzian  in  the  asymptotic 
region.  For  large  v  it  is  given  by 


7  2  ft,, 

A„-—=e 

V  277 


(7b) 


For  the  hyperbolic  secant  puls,  he  transition  art 
plitude  van  Tiles  for  pulse  areas  j.,,  n~.  n  integral, 

for  all  det'  lings.  The  zeros  of  a:l  .  on  the  other 
hand,  occur  for  ft,  rnr  for  zero  detuning,  while 
those  for  large  detuning  are  ft,  2n  1  hose  ot  .... 
go  from  n  ~  at  a  -  0  to  n  ~/2  as  a  -  /. . 

The  existence  of  a  pole  a!  x  i  is  a  sufficient,  but 
not  a  necessary,  condition  that  the  asymptotic 
Fourier  transform  of  a  coupling  pulse  vary  as 
pUo)e~  "  .  For  example,  the  function  (I  -t-.v‘)  '/‘ 
has  an  asymptotic  Fourier  transform  proportional  to 
vl/:e^v.  The  factor  vl/2  precludes  deducing  the 
asymptotic  transition  amplitude  from  the  Roseti- 
Zener  formula.  Similarly,  the  squares  of  the  hyper¬ 
bolic  secant  and  of  the  Lorentzian  each  have  poles 
of  second  order  at  x  =  i  with  the  consequence  that, 
for  both  of  these,  Aa  ~v'e~  '  ,  so  that  while  these 
will  have  asymptotic  transition  amplitudes  that  are 
related  to  each  other,  they  cannot  be  obtained  by 
scaling  from  Eq.  (4).  In  a  future  paper,  we  shall 
show  how  to  calculate  asymptotic  transition  ampli¬ 
tudes  when  the  coupling  pulse  has  second-  and 
higher-order  poles  at  x  =  /.  For  now  ,  we  merely 
present  the  formulas  for  the  transition  amplitudes 
generated  by  the  squares  of  the  hyperbolic  secant 
and  Lorentzian 


Since  the  Rosen-Zener  solution  gives  the  transi¬ 
tion  amplitude  for  all  detunings,  according  to  Eq. 
(4),  as  (recall  that  A  —ft} =Sf) 


we  have  shown  that  the  asymptotic  Fourier 
transforms  of  the  Lorentzian  and  hyperbolic  secant 
are  proportional  for  large  detunings,  the  Lorentzian 
must  induce  a  transition  amplitude  that  obeys  a  for¬ 
mula  similar  to  Eq.  (4).  From  Eqs.  (7),  we  see  that 
to  construct  the  Lorentzian  and  hyperbolic  secant 
Fourier  transforms  so  that  they  are  asymptotically 
identical,  it  is  necessary  to  choose  the  Lorentzian 
pulse  area  fti_  to  be  twice  that  of  ft,,.  Since  f „  —  2 f, 
and  ft,,  —  ft,  /2,  the  asymptotic  transition  amplitude 
for  the  Lorentzian  pulse  may  be  obtained  from  the 
known  result  for  the  hyperbolic  secant  pulse  as 


a,(//2)=-;^e 

c- 


sin 


C 


i  aft 


—  iv  2rrfn(a  )sin/3/y  , 

Xsinh 

C 

_«/L 

1/2 

this  formula  must  be  valid  asymptotically  also.  As 

~ 

X  sinh  C 


.jo 'Pi 

2?r 


iz: 


(9a) 


a -,{L  2)  =  -i2-':e~ 

a  sin 

jc 

aft 

)/: 

c- 

I  1 

2  77 

ii/:  ] 

(9b) 


where  C  =  1  +  -(-•  ■■  -- 1.198.  Equation 

(9a)  can  be  obtained  from  Eq.  (9b)  by  scaling  tech¬ 
niques  derived  in  this  paper.  Equation  (9ai  ts  valid 
only  for  i  ft ,  <  I a  ,  and  Eq.  (9b)  for  (l  <  2<r  . 


a zi  —  —i v  2v  2fi  (a  )sin  -  (8a) 

This  result  has  been  independently  obtained  by  car¬ 
rying  out  an  asymptotic  solution  of  Eqs.  (3).12  One 
can  also  show  that  for  the  pulse 

Ac  -  ft cx  cosechTr.v  , 

the  appropriate  scaling  law  is 

v  2rr  7  .  „  , 

=  i  ^  fftn  )s\n2ftr  .  t Kb) 


III.  SUMMARY  AND  CONCLUSION 

In  this  paper,  we  have  demonstrated  that  pulse 
shapes  Aft)  whose  Fourier  transforms  asymptotical¬ 
ly  approach  the  form  cM  v  le  '  ,  where  6  is  slowlv 
varying,  may  be  categorized  into  families  which 
differ  according  to  the  function  <!>.  Within  each 
family,  the  transition  amplitudes  a:  i  /  1  are  related 
bv  simple  scaling  laws,  so  that  if  one  is  able  to 
derive  an  expression  for  the  transition  amplitude 
generated  by  one  member  of  the  family,  correspond 
mg  formulas  foi  all  other  mcmbcis  of  the  tamilv 
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may  be  written  down  by  inspection. 

A  sufficient  condition  that  the  Fourier  transform 
be  of  the  required  form  is  that  it  be  obtainable  in  the 
asymptotic  region  as  a  contour  integral  evaluated 
from  the  residue  at  a  single  pole  on  the  imaginary 
time  axis.  For  the  case  where  .-11/1  has  simple  poles, 
o  (sc)  may  be  inferred  from  the  solution  of  the 
Rosen-Zener  problem, :  '  known  for  50  years,  by  a 
trivial  scaling  operation. 

Our  results  were  obtained  by  examining  the  struc¬ 
ture  of  the  terms  in  perturbation  expansions  for 
transition  amplitudes.  (We  have  demonstrated  that 
these  sequences  always  converge  in  two-level  prob¬ 
lems  provided  that  the  pulse  areas  are  finite.  Low- 
order  approximations,  however,  are  frequently  not 
useful  for  t  -*  x  even  when  they  are  valid  at  finite 


times. I  With  suitable  of  ratios  of  pulse 

areas,  corresponding  terms  in  the  series  for  different 
members  of  the  same  family  will  be  identical. 

In  a  future  paper,1"  we  shall  present  methods  for 
explicitly  calculating  transition  amplitudes  that  ap¬ 
ply  to  higher-order,  as  well  as  simple  poles  Thus 
we  are  not  restricted  in  practice  to  writing  scaling 
laws  for  pulses  which  may  be  compared  in  the 
asymptotic  region  to  the  hyperbolic  secant. 
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APPENDIX:  CONVERGENCE  OF  PERTURBATION  THEORY  FOR  THE  TRANSITION  AMPLITUDE 


We  demonstrate  here  that  the  perturbation  series  for  a:  converges  for  all  finite  pulse  areas.  The  contribution 
of  order  (2 k  +1)  is 

b\k  +  ]  a  {k  +  ' 

Ik  I 

,,  -  tax  i  ,  w-  -» 

/(.V|)e  dx,  [] 

i  =  - 

Now  assume  that  A  Lx)  is  of  a  single  algebraic  sign.  Without  loss  of  generality  we  may  take  this  to  be  positive. 
We  compare  the  series  with  the  corresponding  expansion  for  a  =  0, 


J  -  X 


/ 


/-I 


..  I  ■  -  I  Jax ,  , 
JtXj  )e  Ul.Xj 


(All 


b'w=  -ip:k+'(  -  1)*  f  /<Jc i  >dx,  FI  f  'X'f(xj,dxj 

j--2 

=  -i(32k  +  ’(  -  1  )*  f  /Lx,)  d.x  |  {J  /  fLXj)  d.Xj  . 


Invoking  the  theorems  on  repeated  integrals  of  the  same  function 

i/3:k  +  1  .  i  /»  *  i  -*  + 1 


b"]  ~  72 k  -t  1 )! 


1-1)*  f  f  Lx  )d.x 


(A2) 

<A2'I 


and  the  terms  are  recognized  as  identical  to  those  for  the  series  i  sm/i.  Now  consider  the  series 


F<P)=  2 


b  in 


J 

(2  k 


:k  t- 1 


J~  f(.x)dx 


,:k  ^  i 


J*  . 

(2A  -r  1  )! 


This  is  evidently  the  series  for  sinh/3,  which  converges  as  long  as  /J  is  finite.  Hence,  the  series  of  Eq.  (A2)  is  ab¬ 
solutely  convergent.  Now 


}'-k  - 1 


„  2k  4-  1  _  X  .  ,  ,i 

f  /(-x,)e  mX'dx  j  n  J"  “X,dxj 

j 


<  p:k+'  fx  !/(*,)■«&,  n’  fM'-,i/<xj)!dxj 

J  X  *  ^  —  CO 


so  that  the  series,  F.q.  (AH  is  also  absolutely  conver¬ 
gent,  and  our  result  is  established. 

We  note  that  the  same  arguments  will  apply  to 
perturbation  scries  at  finite  times,  provided  merely 


that  f  ftx’idx'  -  ftLx)  is  of  one  sign  and  finite. 
If  /  lx  I  changes  sign,  the  results  will  still  be  valid 
provided  the  generalized  area  J  / i.v'i  d\  is 
finite 
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A  simple  case  where  the  convergence  theorem 
does  not  apply  is  the  coupling  function 

A  (x)  =  (const  Htanh7rx/2)/x  , 

since  P  is  logarithmically  divergent.  In  addition, 
since  the  pulse  area  is  proportional  to  the  Fourier 
transform  at  zero  frequency,  the  multiple  integrals 
in  the  frequency  domain  for  the  third-  and  higher- 
order  contributions  to  the  perturbation  series  contain 
regions  where  the  integrands  blow  up,  so  that  the  in¬ 


dividual  terms  beyond  first  order  may  not  even  ex¬ 
ist.  (The  first-order  contribution  will  be  finite,  since 
the  Fourier  transform  for  this  pulse  exists  for  v  0. 
In  this  case,  we  note  that  the  infinite  area  does  not 
imply  a  pulse  of  infinite  energy,  so  that  it  theoreti¬ 
cally  could  exist.  One  evidently  cannot  use  the 
methods  developed  here  to  describe  the  dynamics. 
At  the  very  least,  decay  would  have  to  be  included 
in  the  analysis,  and  a  completely  nonperturbative 
treatment  utilized.) 
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Abstract 

X  r 

With  acci^ite  experimental  values  of  energies,  transition 
moments  and  decay  rates  as  input,  our  calculated  results  show 
that  photon-catalyzed  bound-continuum  processes  can  be  ob¬ 
served  for  the  first  time  using  ^  in  a  cell  or  a  molecular- 

7  2 

beam  experiment  at  laser  intensities  as  low  as  10  W/cm  and 
5  2 

10  W/cm  ,  respectively.  Two  new  features,  post-saturation 
quenching  of  fluorescence  and  huge  (  > 150-fold)  enhancement 
of  fragments,  are  shown  for  the  first  time. 
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Interest  in  laser-induced  multiphoton  bound-continuum  processes 


(e.g.  ionization,  autoionization,  dissociation  and  predissociation) 

1  fkt.  2 

in  gases  and  condensed  phases  continues  to  grow.  In  view  of  the 

tecerdly 

paeefit  high  nonlinearity  (direct  33-photon  ionization). observed  in 

A 

multiphoton  bound- continuum  processes  with  net  absorption,^-  it  is 
surprising  that  there  is  so  far  no  experimental  demonstration  of  the 
lowest-order  (two-photon)  process  in  another  kind  of  multiphoton 
bound-continuum  transitions, 

AB  +  Nhcc  -*•  [ AB '  +  (N±l)hu)]  -*•  A  +  B  +  Nhu  ,  (1) 


which  conserves  the  total  number  N  of  photons  in  a  given  mode  of  the 
electromagnetic  field.  In  reaction  (1) ,  A  and  B  refer  to  neutral  or 
charged  fragments,  and  the  square  brackets  denote  an  intermediate 
state.  An  example  is  given  in  Figure  1.  Since  the  rate  in  the 
reaction  (1)  is  enhanced  by  the  photons  without  their  being  consumed 
in  the  transition  (while  there  are  usually  some  losses  or  gains  due 
to  other  unavoidable  processes  present  in  an  experiment),  this 
multiphoton  process  has  been  called  the  photon-as-catalyst  effect 

3 

(PCE) .  It  can  be  viewed  in  general  as  a  bound-continuum  state  mixing 

due  to  the  external  coherent  .adiation  and  as  a  nonstoichiometric  use 
4-8 

of  laser  photons.  When  a  resonant  intermediate  discrete  level  in 

4  6 

absorptive  multiphoton  processes  (AMP)  is  imbedded  in  a  vibrational  ’ 

5  7  8 

or  ionization  ’  ’  continuum,  PCE  may  occur  simultaneously  (and  some¬ 
times  unsuspectedly )  with  the  AMP. 

3-8 

The  considerable  amount  of  recent  theoretical  works  have 
heightened  the  interest  for  an  observation  of  the  PCE.  To  facilitate 
this  goal,  we  propose  here  two  such  experiments  and  show  their  feasi¬ 
bility  with  detailed  analyses  on  the  bound-continuum  process:  photon- 

« 

catalyzed  predissociation  of  I2  in  selectively  excted  BO+vJ 


states,  as  explained  in  Fig.  1.  The  unpopulated  high-lying  XO  v'J' 

S 

4 

states  are  chosen  as  intermediate  states.  Our  calculation  shows 
for  the  first  time  that  the  PCE  could  cause  a  new  characteristic 
post- saturation  quenching  in  the  BvJ  fluorescence  (a  phenomenon 

new  to  AMP  as  well)  and  a  huge  (  > 150-fold)  enhancement  of  the  atomic 
iodine  fragments.  Our  results  indicate  that  the  PCE  induced  by  1-ys 
laser  pulses  could  be  observed  in  a  cell  experiment  monitoring  the 

7  in 

fluorescence  with  laser  intensity  1  s=  10  W/cm  or  a  molecular-beam 

A 

5  2 

experiment  monitoring  the  fragments  at  I  s  10  W/cm  .  This  range  is 

1112  2 

many  orders  of  magnitude  lower  than  the  requirements  ( e»10  -10  W/cm 

7  8 

calculated  for  other  systems.  ’  In  this  note  we  describe  these  key 
findings . 

Our  calculation  is  based  on  the  published  analytic  results  (Eqs. 
(2.28)  and  (2.30)  of  Ref.  5)  of  a  theory  of  resonant  photon-catalyzed 
bound-continuum  transitions  (via  a  bound  intermediate  state  as  illus¬ 
trated  in  Fig.  1)  regardless  whether  the  continuum  is  vibrational 
or  electronic  and  whether  the  system  is  in  the  gaseous  or  the  condensed 
phases.  The  main  approximation  is  that  the  laser  is  a  square  pulse  of 
constant  intensity  I  and  duration  T.  According  to  Section  III  of  Ref. 
5,  some  favorable  conditions  for  observing  the  PCE  are  (1)  resonant 
intermediate  states;  (2)  large  transition  moments;  (3)  long-lived 
bound  states  (and  long  laser  pulse)  for  long  interaction  time;  and 

rede  5 

(4)  small  competing  spontaneous  bound-free  transition^  .  Indeed  Ij 

9 

possesses  these  advantages:  (1)  the  X-B  transition  frequencies  being 
in  the  tunable  range  of  narrow-band  dye  lasers;  (2)  favorable  elec¬ 
tronic  transition  moments  and  Franck-Condon  factors  and  densities  ; 

(3)  long  lifetimes  (^  lys)  of  BvJ  and  Xv'J'  states^-  and  (4) 

weak  natural  predissociation  (^  10^s  ^). 


These  explain  the  low 


intensity  requirement  and  huge  fragment  enhancement  in  our  results. 


In  contrast,  the  other  calculated  systems  suffer  the  disadvantages 

of  shorter-lived  states,  large  competing  spontaneous  autoionization 

g  7 

rate  and  far  from  resonance  or  small  transition  moments. 

Our  first  set  of  results  is  applicable  to  an  ^  cell  experiment 

in  which  a  resolved  BvJ  -*•  Xv^J^  fluorescence  line  is  monitored.  The 

BvJ  level  can  be  selectively  populated  from  the  ground  state  by  a 

dye-laser  pump  pulse. ^  A  second  time-delayed  dye-laser  pulse  (the 

PCE  laser)  is  incident  on  the  excited  Ij,  and  the  monitored  fluorescence 

line  is  integrated  over  the  PCE-pulse  duration  T.  The  parameters  of 

the  level  system  defined  in  Fig.  1  (chosen  to  match  the  optimum  power 

output  of  laser  dyes  and  for  their  favorable  Franck-Condon  densities) 

are:  XVqJq  =  X’  Wp  =  1?>688.936  cm  1;  BvJ  =  B,  18,  37;  = 

17,475.806  cm-1;  Xv^  *  X,  1,  38;  Xv’J'  =  X,  7,  36  and  the  BvJ-Xv’J’ 

frequency  u)q  »  16,228.577  cm  1.  A  first  step  in  the  experiment  is  to 

tune  the  PCE-laser  frequency  w  into  resonance  with  u)q.  Our  calculated 

results  on  the  integrated  fluorescence  signal  (IFS)  versus  to  show  a 

(power-broadened)  dip  at  u)q.  With  w  locked  at  U)q,  the  next  step  is  to 

increase  the  PCE-laser  intensity  1  and  record  the  IFS  versus  I.  Our 

calculated  values  for  the  ^  vapor  pressure  0.36  Torr  (at  300°K) 

are  plotted  in  Fig.  2.  Considering  first  the  resonance  curves  (labelled 

by  0)  for  10-ns  pulses  and  starting  from  the  low  intensity  region, 

we  note  the  decrease  in  the  IFS  due  to  laser-stimui&ted  B  -*•  X  transfer 

of  population,  and  then  the  coherent  saturation  of  the  B  -  X  transition, 

13 

a  known  phenomenon.  For  1-ys  pulses,  the  resonance  IFS  has  already 

2  2 

reached  saturation  at  I  *  10  W/cm  .  Then  at  still  higher  intensity, 
if  the  PCE  were  to  remain  negligible,  the  IFS  would  stay  constant  (the 


flat  dotted  curves  in  Fig.  2).  In  fact,  however,  our  calculations 
(the  solid  lines)  reveal  a  new  feature:  a  laser-intensity-dependent 
decrease  beyond  the  saturation  region.  This  is  due  to  the  fact  that 
the  laser-stimulated  bound-free  absorption  in  the  second  (and  bottle¬ 
neck)  stage  of  the  PCE  (hence  the  whole  process  itself)  dominates 
over  all  the  other  loss  mechanisms  of  the  bound-state  population.  The 
resulting  inflection  region  between  the  two  quenching  regimes  would 
be  a  readily  recognizable  feature  in  an  experiment.  This  inflection 
persists  for  detunings  up  to  about  0.03  cm  \  but  disappears  for 
A  >  0.3  cm  In  the  latter  case,  the  mark  of  PCE  is  quenching  the 
fluorescence  to  vanishing  values  at  higher  intensity  without  ever  reaching 

saturation,  as  shown  in  Fig.  2.  From  the  A  =  0  curves,  we  see  that  the 

8  2  7 

PCE  could  be  identified  at  I  s  10  Vl/cm  for  10-ns  pulses  and  at  Is=10 

2 

W/cm  for  1-ps  pulses. 

If  the  sharp  resonance  dips  at  characteristic  frequencies  Uq  and 
the  quenching  beyond  saturation  were  observed  in  an  experiment,  we  think 
they  would  constitute  sufficiently  strong  evidence  for  the  PCE.  The 
dependence  on  Wq  means  that  the  level  Xv'J'  must  play  a  critical  role. 

Then  from  energetic  and  symmetry  considerations,  the  only  possibility 
of  a  laser-stimulated  decay  from  the  Xv'J'  level,  and  out  of  the  co¬ 
herently  superposed  B-X  states,  is  a  transition  to  the  lu  or  A  dissocia¬ 
tive  states.  The  dependence  on  the  Xv'J'  levels  would  also  eliminate 
the  possibility  of  the  observed  behavior  being  caused  by  (a)  single¬ 
photon  dissociation,  nonresonant  Raman  or  Raleigh  scattering,  all  from 
the  BvJ  level;  or  (b)  resonant  Raman  or  Raleigh  scattering,  or  resonant 
multiphoton  dissociation  or  ionization  via  some  discrete  levels 
higher  than  BvJ,  since  no  higher  levels  possess  the  identical  spectrum 
as  the  B-X  transitions. 


6 


Our  second  set  of  results  is  applicable  to  a  crossed  laser-molecular- 

beam  experiment  with  the  apparatus  similar  to  those  in  double  absorption 

photofragment  spectroscopy,  capable  of  mass,  translational-energy  and 

14 

angular-distribution  resolution.  The  molecular  beam,  the  nearly  coaxial 

pump-  and  PCE-laser  beams,  and  the  quadrupole  mass  spectrometer  are 

oriented  along  the  x,  y,  z  axes  respectively  (see  inset  in  Fig.  3) .  The 

atomic-iodine  fragments  (as  the  signal)  are  collected  during  the  PCE- 

laser  pulse  and  afterwards.  To  eliminate  the  background  noise  arising 

from  XVqJq  -*•  (lu,A)  photodissociation  fragments  due  to  the  weaker  pump 

pulse,  its  linear  polarization  e’  should  be  oriented  along  the 

14 

detection  z-axis  (under  the  axial-recoil  approximation,  fragment  an¬ 
gular  distribution  “  sin^O ,  with  6  measured  from  £').  For  the  same 
purpose,  the  frequency  oj  of  the  time-delayed  PCE-laser  pulse  is  chosen 
to  be  smaller  than  the  dissociation  energy  DQ(  12,440  cm  1)  of  the 
ground  level.  To  maximize  the  PCE-fragment  signal  the  linear  polariza¬ 
tion  c  of  the  PCE-laser  should  be  oriented  at  43°  from  the  z-axis, 

2  2 

since  the  PCE-fragment  angular  distribution  is  given  by  cos  $  sin  $, 

A 

with  $  being  measured  from  e. 

The  level  system  (see  Fig.  1)  for  the  fragment  calculation  is  given 
by:  XvQJ0  =  X,  0,  54;  u  =  17,280.974  cm"1;  BvJ  =  B,  14,  53;  u>f  = 
16,856.283  cm"1;  Xv^  =  X,  2,  54;  U)Q  =  12,168.112  cm-1  and  Xv'J'  = 

X,  26,  54.  The  BvJ  level  is  chosen  for  its  small  predissociation  rate11 
and  the  Xv'J'  level  for  w<Dq.  The  calculated  total  fragments  collected 
(S)  with  the  PCE-laser  on,  showed  a  resonance  peak  at  u  =  Wq.  Observation 
of  such  peak  would  ensure  the  r  itical  role  of  the  Xv'J'  level  and 
eliminate  the  possibility  of  other  processes,  as  discussed  under  the 
fluorescence  results.  Figure  3  gives  our  results  for  the  fragment 
enhancement  ratio,  R  =  S  (with  PCE-laser  on)/S  (without  PCE-laser).  At 


5  2 

resonance,  R  reaches  the  value  2  at  I  as  low  as  10  W/cm  and  the  maximum 
7  2 

value  152  at  I  —10  W/cm  .  At  low  intensity  where  the  PCE  is  negligible, 
notice  that  R  can  be  (slightly)  less  than  1  owing  to  laser-stimulated 
B-X  population  flopping.  Observation  of  the  resonance  peak  and  of  the 
intensity-dependent  enhancement  in  I-atom  fragments  at  the  expected 
mass,  translational  energy,  and  angular  distribution  would  constitute 
direct  proof  of  the  PCE. 

Based  on  our  earlier  preliminary  results  on  the  fluorescence  system, 

an  experiment  to  observe  the  PCE  was  attempted  by  Brechignac,  Cahuzac 

and  Vetter, 16  using  a  narrow-band  (30  MHz),  single-mode,  YAG-pumped 

dye  laser  (t— 10ns).  Only  the  fluorescence  line  (and  not  the  fragments) 

was  monitored  in  their  cell  experiment.  The  predicted  resonance  dip 

was  observed.  In  the  intensity  study,  they  reported  an  IFS  decrease  of 
8  2 

40%  at  2  x  10  W/cm  ,  of  20%  and  10%  at  lower  (but  unmeasured)  inten- 

7  2 

sities,  and  no  decrease  below  2  x  10  W/cm  .  We  find  that  all  their 

results  can  be  accounted  for  by  the  A  =  0.3  cm  1  solid  curve  in  Fig. 

2.  A  detuning  of  0.3  cm  1  might  be  due  to  a  reported  frequency-drift 

problem  in  their  dye  laser.  If  so,  their  data  point  with  the  highest 
8  2 

intensity  (2  x  10  W/cm  )  has  just  reached  the  region  where  the  case 
with  PCE  begins  to  differ  from  the  case  with  only  stimulated  emission 
(see  Fig.  2).  As  such,  their  results  probably  could  not  be  considered 
as  clear  evidence  nor  as  a  disproof  of  the  existence  of  the  PCE  in  I2 
(BvJ ) . 

The  study  of  Brechignac  £t  a_l.  prompted  us  to  check  the  effects 
of  spatial  and  spectral  averagings  of  the  pulses  and  of  optical  Stark 
shifts.  Of  these,  only  spectral  averaging  (over  a  Gaussian  line 
profile  of  FWHM  0.03  cm-1  resulting  from  overlapping  Doppler-broadened 
hyperfine  lines  of  the  B-X  transitions16)  changes  significantly  the 
approach  to  saturation  for  the  A  =  0  curve,  while  leaving  its 
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PCE-dominant  region  and  the  entire  curves  for  A  =  0.03  and  0.3  cm 

essentially  che  same  as  the  -o^averaged  curves  (see  Fig.  2).  Also,  our 

calculations  can  reproduce  quite  satisfactorily  the  shapes  of  the  three 

experimental  saturation  curves  (each  with  5  data  points)  in  the  ion-dip 

13 

spectroscopy  of  From  these  checks,  we  believe  our  present  results 

are  reliable.  We  hope  these  results  would  stimulate  and  aid  experimental 

efforts  to  observe  this  effect.  We  also  hope  that  this  work  would 

stimulate  studies  of  PCE  in  other  areas  of  physics  and  chemistry,  since 
3  5 

theories  ’  (without  restricting  the  nature  of  the  bound  and  the  continuum 

states)  indicate  that  the  PCE  is  also  applicable  to  other  bound-continuum 

7  8 

transitions:  those  involving  electronic  motion  ’  and  those  occuring 

9  17 

in  liquids  and  solids.  '* 
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Figure  Captions 

Fig.  1.  The  photon-as-catalyst  effect  (the  solid-line  arrows)  on  in  a 

BvJ  level  consists  of  the  stimulated  emission  to  an  unpopulated 

intermediate  level  Xv’J'  and  the  absorption  to  the  vibrational 

continua  of  lu  and  A  states,  to,  to  and  to,  are  the  PCE-laser, 

P  f 

pump-laser  and  fluorescence  frequencies,  respectively. 

Fig.  2.  Integrated  fluorescence  signal  versus  laser  intensity 

for  10-ns  (unlabeled)  and  1-ys  pulses,  plotted  relative  to  the 
natural  fluorescence  amount  (normalized  to  the  same  dashed  curve 
for  both  10-ns  and  1-ys  cases).  Each  pair  consisting  of  a  solid 
curve  (case  with  PCE)  and  of  a  dotted  curve  (case  with  i-r  bound- 
bound  transition  but  without  r-f  bound-free  transition)  is  labeled 
by  the  detuning  in  cm  ^ .  The  pair  labeled  by  0  is  for  zero  de¬ 
tuning  with  spectral  averaging. 

Fig.  3.  Enhancement  ratio  of  fragments  versus  laser  intensity 

for  1-ys  pulses.  Each  curve  Is  labeled  by  the  detuning  in  cm  . 
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INTRODUCTION 

The  study  of  two- level  atoms  coupled  by  external  fields  is 
more  than  50  years  old,  but  continues  to  be  of  interest.1-^  Our 
purpose  here  is  to  deepen  the  A~  eigenvalue  method,  an  approach 
to  this  problem  introduced  recently/-' 

Except  for  a  few  potentials,  closed- form  solutions  are  not 
known. 2»/l  Recently,  Bambini  and  Berman0  (B  and  f,)  found  a  class 
of  coupling  functions  for  which  the  two-state  equations  of  notion 
could  he  solved  analytically.  The  set  includes  the  hyperbolic 
secant  of  Rosen  and  Zener  as  a  special  case.  Apart  from  this, 
all  members  of  the  sat  are  temporally  asymmetric.  For  these,  B 
and  B  shored  that  there  are  no  coupling  strengths  where  the 
transition  probability  vanishes,  excent  on  resonance.0  This 
differs  f  rom  the  hyperbolic,  secant,  where  P  =*  0  for  .any  pulse 
area  A  equal  to  an  integral  multiple  of  m,  regardless  of  It 

is  also  known  that  other  temporally  symmetric  pulses  have  pulse 
areas  for  P  *  0  off  resonance. 

Tliis  suggests  that  it  may  be  true  in  general  that  symmetric 
(asymmetric)  pulses  possess  (lack)  nodes  in  P(A).  The  author 
studied  this  point  by  regarding  the  equations  of  notion  .13  an 
eigenvalue  (EV)  problem  for  By  determining  when  these  EV 

werCgrcal  or  comp lex,  he  was  able  to  generalize  the  B  and  B  re¬ 
sult  ,  finding  that  symmetric  pulses  always  have  nodes,  but 
asymmetric  pulses  do  not,  except  under  over-determined  condi¬ 
tions,  ^ 
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In  the  present  paper,  we  direct  our  attention  to  actually 
determining  EV  and  indicating  how  one  may  perform  accurate 
approximate  calculations  of  them.  We  shall  present  examples 
where  this  knowledge  provides  partial  understanding  of  spectra. 

In  addition,  we  shall  demonstrate  how  to  express  transition  ampli- 
tudes  in  terms  of  the  EV  and  eigenfunctions  and,  in  turn,  will 
exhibit  an  approximation  method  for  P. 

THE  EIGENVALUE  PROBLEM  AND  ITS  VARIATIONAL  APPROXIMATION 

The  time-dependent  Schrodir.gcr  equation  for  amplitudes  a^ 
and  &2  is,  in  RWA  with  detuning  A,  ^ 

i  "  V(t)  eiAt  a 2  *  (la) 

i  a2  "  V(t)  e-iAt  aL  ,  (lb) 

for  real  potentials.  Equivalently 

ax  -  (|  +  iA)  'a  +  V2  a  «  0  ,  (2a) 

•• 

a*2  -  (~  -  iA)  a 2  +  V2  a2  -  0  .  (2b) 

Consider  now  envelope  functions  of  a  single,  algebraic  sign,  and 
define  z  «  /  f(t')  dt'  -  h,  where  f  (t)  «  V(t.)/A,  with 
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Eq.  (3a)  is  subject  to  the  initial  condltJmjr,  a^f-b)  “  G, 
a  '  (-*s)  "  -  iAc’',  <>  real  but  arbitrary.  lor  certain  A  , 

4  o 

(the  EV  An)  02^'^^  also  vanishes.  We  have  previously  shown  that 

A 2  is  real  if  f(t)  is  temporally  symmetric,  while  it  is  ordinarily 

complex  for  asymmetric  potentials  and  A  y  0 . J  The  eigenfunctions 

Corresponding  to  A^  nr{.  a  and  b„  in  the  tv.’o  representations. 

n  -  2n  2n 


We  shall  restrict  the  remainder  of  the  discussion  to  potentials 
that  are  symmetric  in  the  time  (and  in  z). 


Eq.  (4b)  resembles  a  one-dimensional,  time-independent 
Schrodinger  equation  for  a  particle  in  a  complex  "potential. 11  Be¬ 
cause  of  the  non-Hermiticity  of  the  "Hamiltonian",  a  generalized 
form  of  the  usual  eigenfunction  expansion  is  appropriate^. ^  In 
particular,  the  correct  normalization  integral  is  / b  dz,  not 


particular,  the  correct 

/j_Jb  |  r’  dz,  and  matrix  elements  of  operator  0  are  given  by 

0  “*  -n/S,  b  Ob  dz.  Then  A"  is  given  by  (b  normalized  to  1) 
mn  -%  mn  n  b  n 
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If  b  is  exact,  Eq .  (5)  gives  A*"  without  ejjror.  One  may  also  in¬ 
terpret  it  as  a  variational  principal  for  A^.vith  b  now  a  trial 
function,  which,  if  arbitrarily  flexible,  gencratcsnan  Euler- 
Lagrange  equation  that  is  Eq .  (4b).  If  b  is  a  function  of  a 

finite  nv.'ber  of  parameters,  Eq.  (5)  provides  a  variational  ap¬ 
proximation  for  the  eigenvalues  of  Eq.  (4b). 


APPLICATION’  OF  THE  EIGENVALUE  METHOD 


In  this  section,  we  shall  demonstrate  how  the  eigenvalue 
method  can  predict  qualitative  features  of  two-level  spectra.  We 
address  the  question  cf  Che  validity  of  the  Rosen-Zener  conjec¬ 
ture^  fo*  A  small.  For  the  hyperbolic  secant 

P  -  2rr  | V(A)  |  2  sin2  A/A2,  (6) 

where^V  is  the  Fourier  transform  of  V(t)  «  A  f(t).  Boren  and 
Zener  surmised  that  Eq .  (G)  might  be  true  for  all  smooth  V(t), 
replacing  only  the  Fourier  transform  of  the  hyperbolic  secant 
with  ifnat  of  V.  Since  the  result  holds  for  A  =  0.  one  might  ex¬ 
pect  n  region  cf  approximate  validity  for  symmetric  pulses  at 
small  detunings.  It  can  be  established  that;  this  Is  indeed  the 
case,  in  the  sense  that  corrections  are  0(//‘)  if  V  is  differ¬ 
entiable  at  A  «  0.  For  the  Lorentzian,  V  is  not  di.f f ercntiahlc  at. 
A  «*  0.  By  calculating  an  eigenvalue  for  that  potential,  the  it—  Z 
conjecture  will  be  shown  to  be  invalid  for  r;naJ  1  A  In  that  case. 


On  resonance,  the  first  node  in  the  transition  amplitude 

occurs,  for  any  pulse  shape,  at  A  tt.  For  the  purpose  at  hand, 

it  is  sufficient  to  demonstrate  that,  for  small  A,  there  is  a 

linear  correction  in  this  first  eigenvalue  for  the  Lorentzian. 

Off  -resonance,  we  choose  a  trial  eigenfunction  which  reduces  to 

the  known  resonant  b  as  A  0.  This  io 
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where  N  in  a  normalizing  factor,  and  pass  to  the  limit  a  -+■  0. 
This  leads  to 


IT 
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a  correction  to  R-Z  linear  in  A.  A  leading  correction  to  the 
eigenvalue  that  is  linear  in  A  cannot  be  the  signature  of  a  lead¬ 
ing  quadratic  correction  to  the  transition  amplitude.  This  crude,, 
trial  function  gives  a  result  that  is  an  excellent  approximation 
to  the  true  eigenvalue  for  A  =*  0.1. 


AN  EIGENFUNCTION  EXPANSION  FOR  TRANSITION  AMPLITUDES 


Since  the-  b  are  complete,  one  may  use  then  as  a  basis  for 
expressing  unknown  functions.  In  particular,  if  b0  is  the.  solu¬ 


tion  to  Eq.  (4)  for  A‘-  /■  A^ ,  we  r.ny  expand  in  ter:.';  of  the  b 
for  <  z  <  Then  b0(h)  may  be  deduced  f/on  the  expansion  by 


standard  procedures. 
Thus 
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Equations  for  b  ,  b„  are 
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Multiply  Eq .  (8a)  by  b^,  Eq.  (8b)  by  b^,  subtract  and  integrate 
over  all  allowed  z.  This  yields,  using  Green's  thecrem. 


end 
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Since  first  order  theory  is  also  exact  for  all  A  as  z  +  -  ^ 
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A  summation  of  the  form  Z 
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- - -  is  the  Green’s  function 
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For  problems  vhose  exact  solutions  are  not^novn,  G  may  be  npproxi 
mated  by  means  of  a  variational  principle.  * 
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CvliistJnally  aided  radiative  excitation  of  atoms  using  »  laser  detuned  to  the  bine  is  shown  to  lead  to 
heating  of  the  external  degrees  of  freedom.  The  velocity  distribution  of  the  excited  atoms  is  studied  ex¬ 
perimentally  and  computed  with  theoretical  predictions. 
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The.  combined  effect  of  collision^  and  radiative  interac¬ 
tions  upon  an  atom  in  a  vapor  b:..  been  the  subject  of  a 
pre..t  deal  of  theoretical  end  experimental  activity  in  the  past 
several  years.1, 2  A  class  of  such  reactions,  commonly  re¬ 
ferred  to  as  collision.:!!;.'  aided  radi  .’.r  ■_  excitation  (CARE), 
can  be  described  in  a  very  general  v/cy  as 

A  i  -t-  B  i  *f  h  f)  -  A  2  +  B  \  , 


case  will  be  essentially  identical  to  the  Maxwellian  velocity 
distribution  in  the  ground  state;  thi  ,  provides  a  control  mea¬ 
surement  to  which  the  velocity  changes  obtained  in  cod,  ..on 
with  the  heavier  ntr.;  w  arr.O'i  and  xenon,  can  be  com¬ 
pared.  In  the  latter  cases,  the  vdc.cia..:;  of  the  At-  excited 
sodium  atoms  will  be  nontheu  t.d  and  char.»ct''t;/ed  l\  a 
nonequilibrium  dis'tibutton  p(v)  which  has  been  calatlaicd 
using  standard  methods4: 
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power  broadening  (the  natural  width  and  hfs  can  be  neglect¬ 
ed  here).  Study  of  the  dependence  of  the  widths  upon  laser 
power  shows  that  the  power  broadening  is  significant.  How¬ 
ever,  at  a  laser  power  of  50  ntW,  which  was  used  in  all  the 
experimental  runs,  it  is  rather  small  and  is  accounted  for  in 
our  data  processing  at  the  same  time  as  pressure  broaden¬ 
ing,  as  shown  below.  For  the  analysis  of  pressure  broaden¬ 
ing  effects,  we  rely  on  previous  experimental  studies4,9 
which  indicate  that  the  pressure  broadening  coefficients  for 
helium,  argon,  and  xenon  are  the  same  within  15%  for  the 
3/*— 4D  transition.  We  have  thus  assumed  the  same 
Lorentzian  widths  for  the  heavy  rare  gases  as  for  helium  at 
the  same  pressure,  and  we  have  used  the  values  of  the 
Lorentzian  widths  experimentally  determined  by  the  decon¬ 
volution  of  the  curves  obtained  with  helium  to  interpret  the 
absorption  curves  obtained  with  argon  and  xenon.  We  ob¬ 
tained  the  Doppler  contributions  to  the  linov.idths  for  the 
latter  two  passes  by  a  deconvolution  also  assuming  a  Voigt 
profile  (this  assumption  will  be  discussed  below).  Thus  we 
obtained  values  for  the  Doppler  widths  which  are  no  longer 
dependent  on  the  foreign  gas  pressure. 

These  ere  shown  in  Tabic  11.  along  with  the  theoretical 
predictions.  For  both  argon  and  xenon,  it  can  be  seen  that 
the  experimental  values  are  significantly  larger  than  the 
thermal  Doppler  widths.  This  gives  evidence  for  the 
predicted  “heating”  effect.  Nevertheless,  Use  effect  is 
smaller  than  theoretically  calculated.  One  can  account  for 
th'.  if  or..-  considers  the.  processes  that  tbcrrri'dizc  the  -toms 
in  the  3/“  state  Fiist,  once  excited  to  the  .-rate  and  be¬ 
fore  absorbing  a  second  laser  photon,  the  sodium  atoms 
may  undergo  collisions  with  the  perturbers  which  tend  to  re¬ 
turn  the  velocity  distribution  to  the  thermal  one.  V/e  note 
that  the  laser  power  density  is  larger  enough  for  the 
'SP-—4D  tr '.riMtioa  to  be  partly  saturated,  which  reduces 
these  velocity  changes.  Second,  the  rcabsorption  of  reso¬ 
nance  photons,  which  provides  an  alternative  mechanism  to 
excite  atoms  from  (he  ground  state,  has  the  same  effect  of 
washing  out  the  particular  velocity  distribution  of  the  3 /’ 
atoms  by  CARE.  Although  we  work  at  densities  of  sodium 
where  the  optical  depth  for  resonant  light  is  high,  the  transit 
time  of  an  atom  excited  by  CARE  through  the  laser  beam  is 
small,  so  that  rcabsorption  events  take  place  mainly  outside 
of  the  laser  beam,  where  we  do  not  detect  them.  These  cf- 


"Pc.rmaneiH  address:  Laborntoire  dc  Spectroscopic  Hcrtzienne  dc  f 
Ecolc  Nornudc  Supeneure,  l.'ni.crsite  Pierre  cl  Marie  Curie,  4 
Place  Jussieu.  F-75230  Paris.  Ccdex  05  France. 

'J.  L.  Carlstcn  and  A.  SztHe,  Phys.  Rev.  Lett  36,  667  (1976);  J.  L. 
Carlsttn,  A.  Szuke,  end  M.  G.  Rayrner,  Phys.  Rev.  A  j_5,  1029 
(1977). 

JS.  Yeh  an!  P.  R.  Berman,  Phys.  Rev.  A  19,  1105  (1979);  P.  R. 
Berman  end  F.  J.  Robinson,  in  Pl.nton  Assisud  Collisions  and  l{r. 
hted  Tofi.s,  edited  by  N.  K.  Ragmen  end  C.  Guidoni  djarwood- 
Acedernic  Publishers.  Chur,  Switzerland,  1932).  pp.  15-33. 

JS.  Reynaud  and  C.  Cohen-Tannoudji,  J.  Phys.  (Paris)  £3,  1071 
(1932). 

4P.  R.  Berman  anti  S.  Stcnholm,  Opt.  Commun.  24,  155  f  1 978) . 


TABLE  II.  Comparison  of  the  experimentally  obtained  widths 
with  theory  (power  and  pressure  broadening  extracted).  All  values 
are  in  GHz. 


Thermal  Doppler  Theoretical  width  Experimental 
Gas  width  including  heating  Doppler  width 


feels  are  very  difficult  to  estimate  quantitatively.  In  order 
to  study  them  it  would  be  necessary  to  perform  experiments 
at  lower  sodium  and  rarc-gas  pressures,  which  we  were  un¬ 
able  to  do  because  of  signnl-to-noi.se  ratio  considerations. 
The  partial  thcrmalization  in  the  3 P  level  may  also  explain 
why  the  experimental  absorption  line  shape  (big.  4)  for  the 
IP  —  4D  transition  differs  from  the  theoretically  predicted 
one  (Fig.  2),  and  justifies  our  use  of  a  Voigt  profile  in  fit¬ 
ting  the  data. 

In  conclusion,  our  experiment  has  shown  evidence  for 
significant  “heating"  associated  with  the  ccllisionnlly  aided 
radiative  excitation  of  the  3/’ atoms  using  a  laser  detuned  to 
bo  blue.  If  a  !,*•  r,c  enough  fraction  of  the  coins  in  the  User 
beam  undergoes  such  a  transition,  this  effect  might  lend  to 
measurable  Umperaiurc  variations  in  the  wpnr,  providing  a 
new  method  o*'  laser  hearing  or  cooling.  .  i.di  to  p  unt 

out  that  cot.br.a  cannot  he  observed  in  the  ev<v-d  mate,  at 
Ic.-ist  after  one  absorption  process:  If  the  !  ::cr  i  tuned  to 
the  ted  side  of  the  line,  only  fast  atoms  can  he  ext  bed  i  ron: 
the  ground  state,  and  are  Unwed  down,  ivuil'en  in  a  veloci¬ 
ty  distribution  which  can  on  shown  to  be  i.-antic:-.!  to  the 
Maxwellian  one.1''  In  this  case,  she  ground  sMir  is  cooled  j 
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Cooling  of  Vapors  Using  Collislonslly  Aided  Radiative  Excitation 

E.  Glacobino  and  P.R.  Berman 

Physics  Department,  New  York  University,  New  York,  N.Y.  10003 

Collislonslly  aided  radiative  excitation  (CARE)  is  proposed  as  a 

mechanism  for  cooling  an  atomic  vapor.  With  a  CW  laser  power  of 

2 

1.0  W/cm  and  the  resonant  dipole-dipole  interaction  providing  the  colli¬ 
sion  mechanism,  we  estimate  that  temperature  gradients  of  tens  of  degrees 

16  3 

per  cm  can  be  achieved  at  vapor  densities  of  order  10  atoms/cm  . 

Key  words:  cooling,  collisions,  resonant  broadening,  laser  assisted 
collisions. 

1.  Introduction 

Over  the  past  several  years,  there  haB  been  considerable  Interest  In  colll- 
slonal  processes  that  occur  In  the  presence  of  radiation  fields  Buch  as  reactions 
of  the  form 


A2  +  Bx  +  ttfl  ♦  A2  +  Bx 

(A^ ,  atom  A  in  state  1;  B^,  atom  B  in  state  1),  which  are  commonly  referred  to  as 
colllsionally  aided  radiative  excitation  (CARE)  [l].  Most  experimental  studies 
to  date  have  concentrated  on  measurements  of  either  the  CARE  cross  sections  or  the 
frequency  distribution  of  the  reemltted  light.  However,  it  can  easily  be  seen 
that  CARE  also  produces  a  change  in  the  atoms'  translational  energy  and,  as  such 
may  have  potential  use  as  a  method  for  heating  or  cooling  an  atomic  vapor  [2]. 

*Supported  In  part  by  NSF  Grant  I NT  7921530  and  by  the  U.S.  Office  of  Naval 
Research. 
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Let  us  consider  rsdlatlon  of  amplitude  E  and  frequency  ft  acting  on  a  two- 
level  atonic  system  whose  Bohr  frequency  Is  w  *  fi-A  ( | A j  »  Doppler  width).  As 
a  result  of  the  interaction  with  the  field,  the  atom  may  scatter  a  Rayleigh  photon 
or  be  excited  to  its  upper  state.  The  latter  process,  although  possible  without 
collisions  in  strong  fields,  Is  greatly  enhanced  by  collisions  [lj.  Since  the 
final  energy  level  of  the  atom  is  different  from  the  energy  of  the  absorbed  photon, 
the  energy  difference  must  be  compensated  by  a  change  in  the  translational  energy 
of  the  colliding  atoms.  As  pointed  out  previously  [2]  the  net  result  Is  a  cooling 
or  a  heating  of  the  vapor.  Tuning  with  ft  >  u>  produces  heating  and  that  with  ft  <  (d 
produces  cooling. 


2.  Expected  Cooling 


The  magnitude  of  the  effect  depends  strongly  on  the  detuning  A  ■  ft  -  to:  each 

time  an  atom  Interacts  with  a  photon  the  energy  1lA  Is  removed  or  added  to  the  vapor. 

One  may  choose  -hA  to  be  an  appreciable  fraction  of  kT  (at  500°K  kT  corresponds  to 
4  2-1 

a  frequency  of  10  GHz  or  a  wave  number  of  3  x  10  cm  ) ,  so  that  velocity  changes 
are  large  compared  with  those  of  order  1.0  cm/sec  which  occur  as  a  result  of  photon- 
recoil  processes. 

Of  course,  the  process  considered  here  Involves  a  non-resonant  atom-field 

Interaction  and  the  rate  of  excitation  is  reduced  compared  to  a  resonant  excitation 
2  2 

by  a  factor  T  / A  where  T  Is  the  collision  rate.  (For  the  sake  of  simplicity  we 
use  expressions  calculated  using  the  impact  limit  of  line  broadening) .  We  shall 
consider  the  case  of  resonant  broadening;  i.e.  the  resonant  dipole-dipole  inter¬ 
action  between  two  atoms  of  the  species  gives  rise  to  the  rate  T. 

The  rate  of  excitation  to  level  2  from  level  1,  is  given  by 


„  2x2r 


(nx  -  n2) 


(1) 


where  n^  is  the  population  of  state  i,  x  1*  the  Rabl  frequency  defined  by 


(2) 


and  u  is  the  dipole  matrix  element  for  the  transition.  If  n^  »  n2,  the  steady- 
state  upper  population  of  state  2  is 


113 
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(3) 


where  y  is  the  decay  rate  of  the  upper  level.  For  cooling  to  actually  occur,  the 
atom  must  decay  back  to  the  ground  state  via  spontaneous  emission  rather  than  via 
stimulated  emission  (stimulated  emission  would  "heat"  the  atom  if  absorption  cools 
it).  This  means  that  it  is  advisable  to  restrict  the  laser  power  such  that 

<  y  .  (4) 

A* 

beyond  this  limit,  any  additional  laser  power  is  useless  since  there  is  as  much 
stimulated  emission  as  absorption. 

For  an  atomic  density  N  and  cylindrical  interaction  volume  of  length  L  and 
cross-sectional  area  A,  the  power  dH/dt  removed  from  the  sample  is 

NAL,  (5) 

dt  ^2 

which  may  be  rewritten  as 

5?  -  »«•.  <P  <  P„>  <‘> 

where  P  is  the  laser  power  and  P  is  the  maximum  laser  power  consistent  with  con- 

m 

dltlon  (4).  This  result  suggest  a  potentially  large  heating  or  cooling  effect. 
Taking 


•RA/kT  -  0.1  (A/2w  -  lO-*  GHx)  ;  y/2v  -  10  MH*; 


P/P„  *  10 
m 


(P  *  5.0  W/cm2  ;  N  *  1015  cm'3), 


one  finds 


1  dH  -  -1 

NAL  IF  6  ™  8*C 


Under  such  conditions,  the  time  scale  at  which  the  energy  is  removed  or  doubled  is 
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1.0  sec.  In  calculating  dH/dt,  we  assumed  the  ■edlun  to  be  optically  thin. 

The  above  calculation  has  to  be  modified  when  effects  of  radiation  trapping  are 
taken  Into  account.  Equation  (4)  must  be  replaced  by 


(8) 


where  y'  1*  the  effective  upper  state  decay  rate  in  the  presence  of  radiation 
trapping.  Typically,  one  finds  [ 3] 


Y*  *  y/[KL  Or  in  KL)1/2],  (KL  »  1) 


(9) 


where 


IC  -  IW„  (Y/<^) 


(10) 


is  the  absorption  coefficient  for  resonant  radiation  In  a  Doppler  broadened  medium 
(dip  -  Doppler  width).  The  cross  section  oq  Is  given  by 


o 


o 


2U  y? 

eoc*y 


*  A2/2  , 


(ID 


where  y  Is  the  dipole  matrix  element  and  A  the  wavelength  for  the  transition. 
The  rate  of  energy  loss  for  the  sample  Is  now  given  by 


3?  ’  T7y-  f-  <p  <  pi>  <12) 

m 

where  P'  Is  the  maximum  laser  power  consistent  with  condition  (8).  Note  that 
m 

P'  may  be  considerably  less  than  P  .  However,  since  P*/y’  ■  P  /Y»  Eq.  (12)  may 

n  d  m  id 

be  recast  In  the  form 


dH 

dt 


Z_  MAI 

!7y  r NAL 


[P  <  P'  -  p„(y‘/y)] 

n  d 


(13) 


Equation  (13)  Is  Identical  to  the  result  (6)  which  we  found  In  the  absence  of 
radiation  trapping,  except  that  P  Is  limited  to  a  smaller  value. 


Even  though  P  Is  Halted  to  a  value  equal  to  P',  significant  cooling  can 

m 

still  be  achieved.  For  example,  for  the  parameters  chosen  above,  one  can  use 

Eqs.  (9-11)  to  estimate  P'/P  •  y'/y  *  10  6,  assuming  a  cell  length  on  the  order 

■  ■  -7  _1 

of  5cm.  Since  P/P  was  taken  equal  to  10  ,  one  finds  P  s  10  P’  which  does  not 

m  m 

violate  the  condition  P  <  P\  Thus,  the  expected  cooling  is  the  same  as  that 
found  in  Eq.  (7). 

3.  Temperature  Gradient 


We  can  now  estimate  the  temperature  gradient  in  a  cell  created  by  this 
cooling.  Considering  that  the  atoms  are  perfectly  thermalized  on  the  wall  of 
the  cell,  the  amount  of  energy  we  must  remove  per  second,  to  keep  a  temperature 
gradient  dT/dr  between  the  laser  interaction  region  and  walls  is 


I  dH  .  „  £1 

S  dt  dr  » 


(14) 


where  S  is  the  surface  of  the  laser  beam  [S  »  2/nAL  ■  iL;  beam  circumference 
-  £]  and  ic  is  the  thermal  conductivity  of  the  gas  approximately  given  by 


<  *  ku/o 

c 


(15) 


where  k  is  the  Boltzmann  Constant,  u  a  mean  speed,  and  a  a  collision  cross  section. 

c 

Combining  Eqs.  (14)  and  (13)  and  dividing  by  the  wall  temperature  T^  we  obtain 


l_  dT  ^c _ /ISA  \fP  ) 

Ty  dr  “(u/y)lUVVPb/ 


NAL  .  (p  <  p') 

m 


(16) 


We  may  also  note  that  the  limiting  value  P  determined  from  condition  [4]  is  given 
by 


P 

m 


"TcT“ 

o 


A  . 


(17) 


To  optimize  the  temperature  gradient  which  can  be  obtained,  we  proceed  as  follows: 
(1)  We  choose  a  density  N  such  that  the  optical  depth  for  the  laser  radiation  is 


equal  to  unity.  That  is,  w*  set 


V 


-  1  - 


No  Ty 
o  T 


N  -  26  / (o  LTy) 
o 


(18) 


For  L  *  5cm,  and  T  determined  from  theories  of  resonance  broadening  [b]  (l.e. 

r  -  0.023  N  x\),  X  a  600  nm,  and  the  values  of  A  and  y  used  above,  one  finda 
16  3 

N  a  2  x  10  atoms/cm  .  By  choosing  an  optical  depth  of  unity,  ve  ensure  that 

each  photon  "does  its  duty"  in  removing  energy  from  the  sample.  (2)  For  this 

value  of  N,  Eq.  (9)  yields  a  value  y'/y  ■  P'/P_  »  10~2  (3)  Using  Eq.  (17), 

in  m 

we  estimate  the  maximum  allowed  value  P'  for  P  to  be 

IQ 

7s  *  1  X  10"7  Is1  *  0.2  v/cm2  . 

2 

(4)  Since  this  value  of  P'/A  la  readily  obtainable  vlth  a  CW  laser  for  A  <  5  cm 

D 

we  set  P/P  equal  to  its  maximum  allowed  value  (P'/P  )  in  Eq.  (16).  Moreover, 
m  A  mm 

we  take  Noc  “0.1  f/u  to  finally  obtain 


t_  dT 
Tw  dr 


(1  x  10‘7) 


For  A  =  4.0  cm2,  T  -  0.6  x  101®  sec  1,  u2  *  1010  cm/sec,  (liA/kT) 


0.1, 


(19) 


=0.1  (20) 

The  temperature  gradients  predicted  by  Eq.  (19)  (of  order  10  degrees/cm) 
may  be  somewhat  optimistic.  However,  this  order  of  magnitude  calculation  does 
seem  to  indicate  that  significant  cooling  can  be  achieved  using  colllsionally 
aided  radiative  excitation. 

One  of  us  (E.G.)  would  like  to  thank  Dr.  F.  Laloe  for  stimulating  discussions 
on  this  subject. 

*Ue  assume  that  the  excited  state  population  is  close  enough  to  saturation  to 
ensure  that  it  is  resonant  collisions  that  provide  the  major  contribution  to 
the  thermal  conductivity. 
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Over  the  past  several  years  there  has  been  considerable  interest  in  col- 
lisional  processes  that  occur  in  the  presence  of  laser  fields  C 1 3  such  as, 

Aj  +  Bj  +  hfl  ->  A2  +  . 

This  class  of  reactions  is  commonly  referred  to  as  CAFE  (Co.l  1  isional J y  Aided 
Radiative  Excitation).  If  the  laser  is  detuned  from  the  1-2  transition  fre¬ 
quency  by  A,  the  collision  provides  translation  il  kinetic  energy  to  compen¬ 
sate  for  the  mismatch  hA  between  the  field  and  the  atomic  transition,  which 
leads  to  a  corresponding  change  in  the  velocity  <  f  the  atoms  involved  in 
the  process.  This  change  in  translational  energy  is  distributed  between 
the  active  atom  (A)  and  the  perturber  (B)  depending  on  their  relative  masses. 

We  have  studied  the  case  of  the  sodium-rare  gas  system  subjected  to  laser 

award  high,  frequencies, 
tvs  has  been  calculated 
ten  kernel.  The  distribu- 
a  to  the  initial  Maxwellian 
odium* and  takes  practically 
xenon  is  used,  the  sodium'. 


irradiation  detuned  from  the  3S  -  3P  transition 
The  velocity  distribution  of  sodium  excited  3i: 
for  various  rare  gases  using  a  hard  sphere  collar 
tion  obtained  with  helium  perturber s  is  very  clou 
distribution,  since  helium  is  much  lighter  than  r 
all  the  excess  energy.  In  contrast,  when  argon  o 
atoms  arc  predicted  to  undergo  large  velocity  chr.:>g(es ,  arid  the  Doppler  line- 
shape  deviates  from  a  Gaussian  lineshr.pe. 
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Experimentally,  the  vclocit  y  distribution  of  the  3P  atoms  was  monitored 
by  looking  at  the  Doppler  broadened  resonance  associated  with  the  absorption 
on  a  second  transition  3P  *»  AD. 


Wc  observed  a  significant  difference  between  Lbe  Doppler  width  of  the 
3P  -  AD  transition  for  helium,  argon  and  xenon  pm  umbers  (see  Table  I). 

Table  1:  Widths  of  the  absorption  resonance  on  the  probe  3P>AD  in  GHz 
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The  accuracy  is  about  0.1 
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To  precisely  check  by  how  much  the  Doppler  width  w.-ir:  modified  hy  the  heat 
Jng”  effect,  wc  had  to  extract  the  effect  of  |>rc'  : and  power  broadening 
on  the  3)’  :y  AD  probe  transit  ion.  These  were  d(  i  < •, -mined  in  the  ease  of  hcl 
where  there  is  no  heating,  from  the  dcconvolul  ion  <  <  the  expei  imental  Voig 
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profile,  assuming  a  Doppler  width  determined  by  the  Maxwell ian . velocity  dis¬ 
tribution.  We  made  the  assumption  that  the  pressure  and  power  broadening 
effects  leading  to  the  Lorentzian  line  width  were  the  same  for  the  some  pres¬ 
sure  of  helium,  argon  and  xenon  perturbers;  this  is  justified  by  the  fact 
that  the  pressure  broadening  coefficients  on  the  31’  ■*  4D  transition  have 
been  shown  to  be  very  close  to  each  other  for  these  three  rare  gases  [23. 

Then  a  coarse  deconvolution,  assuming  a  Caussian  shape  for  the  Doppler 
broadened  part  of  the  line,  gave  us  the  FWHM  of  this  Doppler  broadened  part, 
to  be  compared  with  the  predicted  theoretical  value  (Table  II)  and  wiLh  the 
width  given  by  the  Maxwellian  distribution. 

Table  II:  Comparison  of  the  experimentally  obtained  width  with  theory 
(power  and  pressure  broadening  extracted) . 


10  There  is  an  obvious  heating  effect  of  the  3P  atoms,  although  it  is  less 

than  theoretically  predicted.  We  interpret  this  discrepancy  as  a  result  of  '  ~ 
re-absorption  of- resonance  photons  which  excite  atoms  from  the  ground  state 

13 — without  changing  their  velocities,  and  of  the  velocity  changing  collisions  - 1 

undergone  by  the  3P  atoms  after  it  has  been  excited.  Studies  at  lower  sodium 

12  _  pressures  are  in  progress  to  check  this  point.  12 

13  In  conclusion,  we  have  demonstrated  velocity  changes  subsequent  to  CARR.  .. 

-  If  the  efficiency  of  the  process  is  high  enough,  that  is,  if  we  are  able  to 

accelerate  or  slow  down  enough  atoms,  this  should  lead  to  a  macroscopic  local 
li —  heating  or  cooling  of  the  vapor,  as  predicted  (33.  '  .  —I'l 

15  This  work  was  supported  by  the  U.S.  Office  of  Naval  Research,  NSF  Grant  - 
INT  7921530,  and  NSF  Grant  PHY  8204402-01. 
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